Momentum and heat transfer from an asymmetrically confined 
circular cylinder in a plane channel 


A thesis submitted 
in partial fulfillment of the requirements 
for the degree of 
Master of Technology 


by 

Srinivas Mettu 



to the 

Department of Chemical Engineering 
INDIAN INSTITUTE OF TECHNOLOGY 
KANPUR 

May, 2005 




1 9 M 

'TRcfVir 


CERTIFICATE 


It is certified that the work contained in the thesis entitled '‘‘'Momentum and heat 
transfer from an asymmetrically confined circular cylinder in a plane channel” by 
Srinivas Mettu, has been carried out under my supervision and that this work has not been 
submitted elsewhere for a degree. 


May, 2005 


Dr. Nishith Verma 

Associate Professor 

Department of Chemical Engineering 

Indian Institute of Technology 

Kanpur 



ACKNOWLEDGEMENT 


At the very beginning I would like to thank the authority of IIT Kanpur for providing 
high-class lacilities and the ambience necessary to carry out research work. During the 
research work, my thesis guide Dr. N. Verma has helped me immensely with his ingenious 
ideas and valuable guidance. I eun highly grateful to him for motivating and encouraging me 
during all the stages of my thesis work. By working under his aegis I have come to know the 
intricate details of research work. It was my privilege to have worked under him. I am also 
thankful to Pro!'. R. P. Chhabra for his kind co-operation and helpful suggestions. 

I would like to specially acknowledge the support regarding CFD provided by 
Ramprakah, Nitin, Dhiman, Sunil Dhole and Nanda Kishore. I shall always cherish the sweet 
incinorics attached with my lab mates Sudhakar, Vivek, Ritesh, Kaushal, Nilesh, Shiva 
Kumar, Rupesh and Abhinandan. They have always helped to maintain a friendly atmosphere 
in our lab. More than this, they have always been a constant source of inspiration. 

The time spent with Anji, Babai, Ravi, Ramana, Simha, Mallesh, Dinesh, Sameer, 
Akhi, Parimal, Satyam, Pramod, tmd Nageswara Rao has been a major recreation for me. The 
memories associated with them will surely be cherished forever. They were always there to 
give me moral support in every situation in my life here. 

At the end 1 would like to convey my regards to my parents and all other family 
members for their endless love, encouragement and endurance during my stay at IIT Kanpur. 
Their silent presence behind my all achievements carmot be properly expressed. 

Finally, I am grateful to the almighty for what I am today. 


May, 2005 


Srinivas Mettu 



ABSTRACT 


I he flow past a circular cylinder is a fundamentally important problem as it is a 
good approximation to a large number of problems encountered in engineering 
applications like heat exchanger tubes, instmmentation technology in Hot-wire 
anemometry, flow past dividers in polymer processing, piping installations. The 
fundamental nature of the flow has attracted many investigators who have carried out 
experimental and numerical investigations to characterize the flow using global 
parameters like; wake length, separation angle, drag coefficient, lift coefficient, Strouhal 
number and Nusselt number. Despite significant number of studies carried out over the 
years on the flow past a circular cylinder in unconfined domain, only few studies are 
available on the flow past a circular cylinder confined in a plane channel that is rather 
more appropriate in many practical situations. In addition, the transition of flow past an 
asymmetrically confined circular cylinder in a plane channel from steady to unsteady 
flow has not been characterized completely compared to the transition of the flow past a 
cylinder in unbounded domain. 

In this study, the unsteady flow of an incompressible Newtonian fluid past a 
circular cylinder in unbounded and bounded domains was simulated using a 
commercially available flow modeling software, FLUENT (version 6.1), which utilizes a 
finite volume based numerical method to solve the equations of continuity, momentum 
and energy. In particular, three cases were studied; (i) unsteady momentum transfer from 
a circular cylinder in unbounded domain (ii) unsteady momentum and heat transfer from 
an asymmetrically confined circular cylinder in a plane channel, and (iii) three- 
dimensional transition in the wake of a circular cylinder. In the case of unsteady 
momentum transfer from a cylinder in unbounded domain, the transition of flow from 
steady to unsteady is studied with the aid of various momentum upwinding schemes. The 
important engineering parameters like total drag coefficient, Strouhal number were 
computed using various upwinding schemes in the range of Re from 50 to 1000 and 
compared with the published data wherever possible. In the case of unsteady momentum 
and heat transfer from an asymmetrically confined cylinder in a plane channel, the 
transition of the flow from steady to unsteady is characterized as a function of blockage 



ratio (/?, ratio of the diameter of the cylinder to the channel width) and gap ratio {y). The 
unsteady numerical simulations were carried out for the range of Re between 10 and 1000 
and blockage ratio (fi) from 0.1 to 0.4 using the constant wall temperature boundary 
condition on the surface of the cylinder. The fluid considered was air with negligible 
property variation. The effect of blockage ratio (p) and gap ratio {y) on the time averaged 
drag coefficient, lift coefficient, Strouhal number, and Nusselt number was studied. In the 
case of three-dimensional flow past a circular cylinder, the transition of the flow from 2- 
D to 3-D flow was studied. 

The critical Re at which the flow transition takes place from steady to unsteady 
was found to increase with increase in blockage ratio (P) for a symmetrically confined 
cylinder in a plane channel. For fixed blockage ratio, the critical Re was found to increase 
with decrease (cylinder moves closer to the wall) in the gap ratio. For fixed Re and 
blockage ratio, the time averaged drag coefficient and Strouhal number were found to 
increase with decrease in gap ratio. The increase in the drag coefficient was found to be 
more pronounced at low Re. The amplitude of the oscillations in the lift coefficient was 
found to increase in the negative direction as the cylinder moved closer to one of the 
confining walls. The reverse trend was observed when the cylinder was present 
sufficiently close to the cylinder in which case the amplitude of the oscillations in the lift 
coefficient increased in positive direction. The effect of increase in blockage ratio on the 
surface averaged Nusselt number was found to be negligible for blockage ratio below 0.3. 
The effect of decrease in gap ratio on the surface averaged Nusselt number was aisio 
found to be negligible. However, some distortion was found in the local Nusselt number 
distribution on the surface of the cylinder when the cylinder was present close to one of 
the confining walls. The transition of the flow past a circular cylinder from 2-D to 3-D 
was characterized using Strouhal-Reynolds number relation and normalized velocity 
profiles along the span of the cylinder. 
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CHAPTER 1 
INTRODUCTION 

1.1 Introduction 

The flow past a circular cylinder is a classical problem in fluid mechanics as it is a 
good approximation to a large number of problems encountered in engineering 
applications like heat exchanger tubes, instrumentation technology in Hot-wire 
anemometry, flow past dividers in polymer processing, piping installations, electric 
transmission lines, suspension bridges, and offshore cylindrical drilling rigs. The 
governing equations of motion, namely Navier-Stokes equations are highly nonlinear in 
nature and the exact analytical solutions are still not available except few approximate 
analytical solutions (Stokes, 1851) for flow past a sphere at low Reynolds numbers, 
rherefore, it has been subjected to numerous experimental as well as numerical 
investigations [1-55] to understand the dynamics of the flow by calculating the global 
parameters that characterize the flow, such as wake length, separation angle, drag and lift 
coefficients, Strouhal number (St), and Nusselt number (Nu^). 

The flow past a circular cylinder may be divided into different regimes depending 
on Re. For i?e ~ 1, the flow is fore and oft symmetric, divides at front stagnation point 
and reunites smoothly behind the cylinder. As Re increases above 5, the flow separation 
takes places on the downstream side with formation of two symmetric vortices, which are 
stable and remain attached to the cylinder. In this regime, symmetric vortices behind the 
cylinder elongate farther as Re increases. Around Re = 47, a new flow pattern develops in 
which the oscillation of vortices gives rise to the well-known Von karman vortex street. 
The discrete vortices are regularly placed behind the cylinder, with the non-dimensional 
number St characterizing the frequency of oscillation of the vortices. In this regime the 
frequency of oscillation (St) increases as Re increases. This regime continues up to Re = 
1 80 at which complex three-dimensional flow structure develops [51]. 

A large amount of wealthy information [1-55] on the flow past a cylinder has 
been accumulated over the years because of the ideal nature of the flow, which can be 
applied to many industrially important applications with little modifications. The above- 
mentioned flow regimes that exist for the flow past a cylinder in unbounded domain may 
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change considerably [23] for the flow past a cylinder confined symmetrically in a plane 
channel. From the literature survey it is evident that the flow past a cylinder confined in a 
plane channel has received relatively much lesser attention [55,23,39], Chen et al. [23] 
have studied the transition of the flow past a circular cylinder confined symmetrically in a 
plane channel from steady to unsteady using finite element numerical simulations. It was 
concluded from their study that the flow dynamics differ greatly when the cylinder is 
confined in a plane channel because of the presence the stationary walls and the stability 
of the flow is enhanced as the walls approach closer to the cylinder. The critical Re at 
which the flow transition takes place from steady to unsteady was found to increase with 
increasing the blockage ratio (/?). 

rhe flow past a cylinder confined in a plane channel is also affected considerably 
depending on the situation whether the cylinder is confined symmetrically in the plane 
channel or confined asymmetrically [24]. In the case of the cylinder confined 
asymmetrically, the flow dynamics changes significantly because of the presence of the 
wall close to the cylinder. To ascertain the effect of asymmetry on the transition of the 
flow, Zovatto and Pedri 22 etti [24] have numerically studied the loss of stability of the 
flow past a circular cylinder confined asymmetrically in a rectangular plane channel for a 
blockage ratio of 0.2 using a finite element numerical method based on vorticity- 
streamfunction formulations. It was concluded from their study that the stability is 
enhanced as the cylinder approaches one of the confining walls. The critical Re at which 
the flow transition takes place was found to increase as the cylinder approaches one of 
the confining walls. The main limitation of their study was that the definition of the 
asymmetry parameter (j, gap ratio) was such that it was applicable to only one blockage 
ratio, which does not allow the comparison of the results with the results for the other 
blockage ratios. In our study we have defined an asymmetry parameter, gap ratio {y), 
such that it takes a value of 1 when the cylinder is confined symmetrically in the plane 
channel, whereas it takes a value of zero when the cylinder touches one of the confining 
walls irrespective of the blockage ratio (fi). We have studied the transition of the flow 
past a cylinder confined asymmetrically in a plane channel over a wide range of blockage 
ratio (/3). More importantly, we have also studied the effect of the asymmetric 
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confinement on the heat transfer characteristics from the cylinder that is more relevant to 
practical applications. 

1.2 Objectives 

The objectives of this thesis work are as follows: 

1 ) Comparison of various momentum upwinding schemes with respect to the transition of 
the flow past a circular cylinder in unbounded domain from steady to unsteady and the 
comparison of global parameters like drag coefficient and Strouhal Number obtained 
using these upwinding schemes. 

2) Characterization of the transition of the flow past a cylinder confined in a plane 
channel from steady to unsteady as a function of asymmetry {y, gap ratio) and blockage 
ratio iP). 

3) Study of the effect of asymmetry and blockage ratio on drag coefficient, lift 
coefficient, Strouhal number, and Nusselt number over the range of Re from 1 0 to 1 000. 

4) Study of the transition of the flow past a circular cylinder in unbounded domain from 
2-D to 3-D. 
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CHAPTER! 

LITERATURE REVIEW 

2.1.Steady flow past a circular cylinder: 

1 he first attempt to solve the Navier-Stokes equations for the problem of steady, 
incompressible, Newtonian flow past a circular cylinder in unconfined domain dates back 
to Ihom [1], who solved the problem by integrating the equations of motions under 
steady state conditions, for Re =10 and 20. Subsequent numerical investigations were 
carried out by Kawaguti [2] and Apelt [3] to solve the same problem for Re = 40. The 
time dependent simulations were first carried out by Payne [4] for Re = 40, 100 and the 
drag values, stream function, iso-vorticity plots were presented. But this study had many 
limitations. For example, a fixed time step was used in the computations and numerical 
computations were carried out for a shorter time at which standing vortex was not fully 
developed. More important, the grid sensitivity analysis was not carried out. In his time 
dependent numerical simulations, Ingham [5] varied time step as well as mesh sizes to 
determine the optimum values for Re =100, Finer mesh size used by Ingham [5] gave 
better results compared to Payne’s [4]. It was concluded that the mesh size used by Payne 
[4] was coarse; hence the steady state values were not obtained. 

Kawaguti and Jain [6] carried out a systematic numerical study of the flow of an 
incompressible, Newtonian fluid past a circular cylinder in unconfined domain for Re = 1 
to 100. The steady state solutions for Re = 10, 20, 30, 40, 50 were obtained as the limit of 
unsteady solutions. The radius of the outer circular boundary was fixed at 1 12 radii of the 
cylinder to have insignificant eflfect of boundaries on solutions. The streamlines and iso- 
vorticity lines were presented for each Re. The length of standing vortex, separation 
angle, total drag, surface pressure and vorticity distribution were plotted as functions of 
time. However, the solutions for Re = 60 and 100 did not reach steady state and it was 
evident from the plot of the length of standing vortex against time for Re = 60 and 100 
that the shape didn’t correspond to steady state flow. It was also opined that the 
appearance of the region of positive (negative) vorticity just outside the wake could have 
been caused either by the coarse mesh used by them or by the occurrence of instability 
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resulting in von Karman vortex street, which appears in the flows at Re > 40. But the 
main limitation of their study was that the grid independence test was not carried out. 

Hamielec and Rail [7] have extended the study of Kawaguti and Jain [6] up to Re 
= 500 with a finer mesh size in radial and angular directions. The radius of the outer 
circular boundary was also varied to study the effect on drag values. The surface vorticity 
distribution, pressure distribution, streamlines, iso-vorticity lines, drag values and length 
of standing vortex were presented for various Re up to 500. It was concluded that the 
length of standing vortex predicted by Kawaguti and Jain was too high owing to the 
coarse mesh size used by them. However, the numerical results for Re = 500 were not 
reliable as the steady state solution was not obtained. In agreement to the observations of 
Kawaguti and Jain (6), fhoman and Szewczyk [8], the authors have noted that at such a 
high Ri\ one experimentally observes vortex shedding and von Karman vortex street, 
which are the characteristics of unsteadiness in the flow. In theoretical studies, vortex 
shedding can only be observed by introducing a small asymmetric disturbance to the flow 
by various techniques such as, introducing artificial vorticity into the flow field for a 
small period of time, rotating the cylinder for a short period of time, and making the 
location of the cylinder slightly eccentric between the two planes. 

Takami and Keller have [9] have carried out numerical simulation of the steady, 
two-dimensional flow past a circular cylinder in infinite domain for Re = \ to 60.The 
radius of the outer circular boundary was varied. The boundary conditions at the outer 
cylinder were applied by matching Imai’s [10] asymptotic solution. Streamlines, vorticity 
plots, drag values and length of standing vortex were presented up to = 60. The linear 
dependence of length of standing vortex was observed up to Re = 60. 

Dennis and Chang [11] extended the study of Takami and Keller [8] to Re =100. 
T he main objective of their study was to produce reliable data on the steady flow past a 
circular cylinder with respect to the growth of standing vortex and to gain an information 
on the nature of the theoretical flow limit as Re-^co. Finite-difference solutions of 
equations of motion were obtained for Re = 5 to 100 for half domain with symmetric 
boundary condition. The length of the standing vortex followed linear dependence on Re 
with the flow separation taking place at a critical i?e of 7. Streamlines, wake length, drag, 
angle of separation, surface pressure distribution, vorticity distribution were presented for 
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Re 5, 7, 10, 20, 40, 70, 100. Although variation in these parameters with Re was found 
to be consistent in their study, the nature of these parameters could not be established 
with certainty as i?e->oo. 

bomberg [12] has numerically studied the steady viscous flow past a circular 
cylinder in unconfmed domain for Re from 20 to 300. He has proposed a new technique 
for the boundary condition at a large distance from the cylinder and an iterative method 
to eliminate the numerical difficulties encountered in the simulation of flows at Re >100. 
'fhe streamlines, surface pressure distribution, vorticity distribution, separation point, 
wake length and drag values were presented. The shortening of the wake region with 
vorticity convected into interior was observed at high Re >260. Far away from the 
cylinder, a simple boundary condition of normal derivatives equated to zero for Re< 40 
and a mixed boundaiy condition for higher Re gave better results. 

D’Alessio and Dennis [13] have proposed a vorticity model for viscous flow past 
a circular cylinder that is consistent with the boundary layer theory for large Re and with 
the asymptotic solution at a large distance from the cylinder for low Re. Drag values were 
presented for Re - 5, 20, 40 and 100 and were in good agreement with the results of 
Dennis and Chang [11] and Fomberg [12]. 

Coutanceau and Bourd [14] have carried out experimental investigation to 
determine the main feature in the wake of a circular cylinder in uniform cross flow in the 
range of Re from 5 to 40 for various blockage ratios Xfratio of the diameter of the 
cylinder to the channel) by flow visualization technique. Velocity distribution in the wake 
of the cylinder, length of the standing vortex, separation angle were plotted against Re for 
various blockage ratios (X). The length of standing vortex and separation angle were 
found to decrease with increase in blockage ratio. The upper limit of Re at which twin 
vortices are stable and adhere to the cylinder was determined to be 39.5. 

Recently, Chakraborty et al. [55] have extensively studied the effect of walls on 
steady, two-dimensional flow past a circular cylinder confined symmetrically in a plane 
channel by numerical simulations for Re between 0.1 and 200 and for blockage ratios (A,) 
between 0.05 and 0.65. Total drag, pressure drag, friction drag, ratio of pressure drag to 
friction drag, the length of recirculation zone and separation angle were plotted against 
Re for various blockage, ratios. For a fixed value of Re, total drag was found to increase 
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with the blockage ratio and for a fixed blockage ratio, drag decreased with increase in Re. 
Piessnre drag and friction drag followed the trend of total drag where as their ratio 
increased with increase in Re for a fixed blockage ratio. The length of recirculation zone 
and separation angle found to decrease with increase in blockage ratio. 

2.2. Unsteady flow past a circular cylinder: 

I'he flow past a circular cylinder may be divided into different regimes depending 
on Re. For Re ~ 1 , the flow is fore and oft symmetric, divides at front stagnation point 
and reunites smoothly behind the cylinder. As Re increases above 5, the flow separation 
takes places on the dovmstream side with formation of two symmetric vortices, which are 
stable and remain attached to the cylinder. In this regime, S 5 anmetric vortices behind the 
cylinder elongate farther as Re increases. Around Re = 47, a new flow pattern develops in 
which the o.scillation of vortices gives rise to the well-known Von karman vortex street. 
The discrete vortices are regularly placed behind the cylinder, with the non-dimensional 
number St characterizing the frequency of oscillation of the vortices. In this regime the 
frequency of oscillation {St) increases as Re increases. This regime continues up to Re = 
1 80 at which complex three-dimensional flow structure develops [51]. 

Many numerical investigations have been carried out to simulate the unsteady 
flow past a circular cylinder. The time dependent numerical simulation of the Navier- 
Stokes equations for the unsteady flow past a circular cylinder in cross flow was first 
carried out by Payne [4] for Re = 40 and 100, followed by Ingham [5] for Re =100. The 
evolution of flow properties like drag and length of the recirculation zone with time was 
presented. 

Collins and Dennis [15] have numerically studied the flow past an impulsively 
started circular cylinder by integrating the time dependent Navier-Stokes equations for Re 
= 40, 100, 200 and 500. The development of various flow properties, namely the length 
of the recirculation zone, the surface vorticty distribution and separation angle with the 
time was presented. For Re < 100, the integration was carried out until the solution 
reached steady state. For Re > 100, the integration was terminated well before the 
solution reached steady state due to the convergence failure. The time at which the 
integration procedure failed to converge decreased with increase in Re. This is consistent 
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witli the fact tliat as Re increases the transition from steady flow to unsteady flow takes 
place at much earlier times. 

1 atel [161 has presented the semi-analytical solutions for the flow past an 
impulsively started circular cylinder for Re =60, 100, 200, 500, 550 and 600. A set of 
coupled parabolic differential equations was obtained by substituting the stream function 
and vorticity expanded in the finite Fourier series in the Navier-Stokes equations. These 
equations were then solved numerically using the finite difference method. The growth of 
the recirculation zone, evolution of drag and surface vorticity distribution was presented 
for each Re. l or the range Re investigated in this study, calculations never reached steady 
state and the length of the wake continued to grow with the appearance of the secondaiy 
vortex on the surface of the cylinder for Re = 500, 550, and 600. 

The main limitation of all the above-mentioned numerical simulations was that it 
was assumed that the flow reached steady state after a long time or the solutions were 
presented for short times at which flow remained steady. But in the flow at Re>40, one 
experimentally observes Von karman vortex street which is a characteristic of 
unsteadiness of the flow. None of the above mentioned authors simulated the vortex 
shedding from the cylinder for Re > 40. 

Braza et ah [17] have solved the time dependent Navier-Stokes equations for the 
flow past a circular cylinder in unconfined domain by a finite volume second order 
accurate scheme and an alternating direction implicit procedure for Re = 100, 200 and 
1000. The radius of the outer circular boundary was taken sufficiently large to simulate 
unbounded flow. They have shown that the flow reaches steady state symmetric pattern 
afler a shorter or longer establishment period even at Re =1000, also observed by 
Martinez. (1979) for Re =100 and 200, attributed to the fact that in case of numerical 
sinmlations the destabilizing factors like non-uniform inlet conditions, irregularity of the 
boundary conditions, vibrations and surface roughness that are present in experimental 
conditions are. absent as the boundary conditions and initial condition are symmetric 
which lead to the steady symmetric solution of the Navier-Stokes equations. That is why 
the authors have used a perturbation method to generate vortex shedding from cylinder, 
which consists of the rotation of the circular cylinder in clockwise direction followed by 
rotation in anticlockwise direction for about a short period of time and then keeping the 
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cylindti stJtiondry for the lest of the calculations. The global characteristics of the flow 
such as the mean drag, evolution of the lift coefficient with time, and mean St were 
piescnted. I he oscillation ol the pressure and viscous drag was found to be in phase with 
the oscillation of the total drag and the mean drag for unsteady flow was higher than that 
for the steady flow. 1 he frequency of the oscillation of stagnation, separation and 
reattachment angle was in phase with the lift coefficient oscillations while the pressure 
coefficient was found to oscillate with twice the frequency of lift coefficient oscillations 
due to the simultaneous contribution of both the alternating vortices. It was shown from 
their study that the secondary eddies appear in the vicinity of each primary eddy above 
Re ■ 200 and have significance presence in the near wake at Re = 1000. 

Jackson [ 18| has studied the onset of periodic vortex shedding from bluff bodies 
of various slrapcs using finite element numerical simulations for two-dimensional laminar 
flow. Time dependent numerical simulations were not performed but the eigenvalue of 
the Jacobian matrix that characterizes the transition from steady to periodic flow was 
calculated by solving the extended systems of equations for locating the symmetry 
breaking Uopf bifurcation, which is the characteristic of the onset of vortex shedding. 
The solutions of these equations gave St and the critical Re at which the flow became 
periodic. The critical Re for the flow past a circular cylinder was determined to be ~ 46 
and it was found to be in good agreement with the experimental studies [19,20] as well as 
with the time dependent numerical simulations [8,21]. The critical Re for the flow past a 
flat plate, ellipse and prism shaped bodies was also presented as functions of the angle of 
orientation with flow. The critical Re was found to increase with increase in the angle 
between the minor axis and flow for an ellipse and with increase in the angle of 
orientation of the flat plate with flow. 

Franke et al. [22] have solved the two-dimensional, unsteady Navier-Stokes 
equations for the flow past a circular as well as square cylinder using a third order 
accurate discretization scheme QUICK for the range of Re < 5000 for the circular 
cylinder case and Re < 300 for the square cylinder case. The main goal of their study was 
to examine the development of the vortex shedding from the cylinder as the flow velocity 
was increased from zero to free stream value. In their study flow velocity was increased 
from zero to the free stream value following a sine function. No disturbance was 
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introduced to obtain vortex shedding. The vortex shedding was automatically triggered 
by lound-olf enois and numerical diffusion. The drag values were in good agreement 
with the calculations of Braxa et al [1 7] for Re up to 300. The deviation at higher Re was 
attributed to the beginning of the stochastic turbulent fluctuations. 

Chen et al. [23] have studied the stability of the flow past a circular cylinder 
confined symmetrically in a plane channel using numerical simulations. Steady flow past 
the circular cylinder was perturbed slightly by the rotation of the circular cylinder for a 
short period of time. 1 he time dependent motion following the rotation was examined to 
determine the critical Re at which the perturbation amplified leading to unsteady flow, 
'fhe critical Rc based on the diameter of the cylinder for a blockage ratio (d/H) of 0.2 was 
found to be ()9. Numerical bifurcation studies, similar to those carried out by Jackson [], 
were carried to determine the critical Re for various blockage ratios ranging from 0.1 to 

0. 7 by the calculation of the eigenvalues of the Jacobian matrix of the flow equation. It 
was showed that stability is lost through a symmetry breaking Hopf bifurcation. The 
critical Re was found to increase with increase in blockage up to 0.5 and decrease with 
the blockage ratio from 0.5 to 0.7. 

Zovatto and Pedrizzetti [24] have numerically studied the loss of stability of the 
flow past a circular cylinder confined asymmetrically in a rectangular plane channel for a 
blockage ratio of 0.2 using a finite element numerical method based on vorticity- 
streamfunction formulations. Flow was assumed to start impulsively and at the inlet a 
parabolic velocity profile was specified. The critical Re at which the flow became 
periodic was found to increase as the cylinder approached one of the two confining walls 

1. c. the vortex shedding from the cylinder was suppressed when the cylinder was closer to 
one of the walls. When the cylinder was placed symmetrically between the plane walls 
the wake resembled the unbounded flow with the symmetric wake. However, when the 
cylinder was closer to one of the walls the wake was asymmetric with a reduced wake 
vorticity on the side of the wall and elongated wake vorticity on the other side. The 
cylinder behaved like a surface mounted obstacle when it was very close to the wall. For 
a fixed Re, the drag coefficient based on the diameter of the cylinder and the non- 
dimensional time period of the oscillations was found to increase with the decrease in the 
gap between the cylinder and the wall. It was concluded that the transition from steady to 
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unsteady flow is delayed as the cylinder approaches one of the confining walls resulting 
in the stability of the flow. 

2.3 Laminar forced convection from a cylinder in cross flow 

The heat transfer from a cylinder is of practical significance in Hot-wire 
anemometry which is a measuring technique widely used in experimental fluid 
mechanics to obtain information on velocity field. The correlation between heat loss from 
the wire probe and the velocity of the fluid helps in determining the flow field. It is an 
indirect measurement technique, which requires calibration against the known flow field. 

Many experimental and numerical investigations have been carried out on 
momentum and heat transfer from a circular cylinder with respect to Hot-wire 
anemometry. Dennis et al. [25] have carried out numerical study for the heat transfer 
from a circular cylinder at constant temperature in a steady stream of viscous, 
incompressible fluid in unconfmed domain in the range of Re from 0.01 to 40 with the 

assumption of constant fluid properties. The authors have determined the variation of Nu 
with Re and compared the values with the experimental results of Collis and Williams 
[26]. The numerical results were in good agreement with the experimental results. They 
have also detennined the variation of local Nu over the surface of the cylinder for the 
range of Re from I to 40. 

To resolve the discrepancy between the data of hotwire anemometers and the 
empirical laws, Hua et al. [27] have theoretically analyzed the heat transfer from a 
constant temperature circular cylinder in cross flow. The authors assumed two regions: a 
stationary fluid region close to the cylinder and potential flow region far away from the 
cylinder. 'Fhe authors theoretically calculated the heat transfer coefficient for airflow over 
the cylinder in the range of Re from 0.01 to 9x10^ and compared the values with the 
experimental data of Collis and Williams [26] and Cole and Roshko [28]. The results 
were in good agreement with the experimental data only in the range of 0,1< i?e < 3x lO'^. 
The deviation at higher Re was attributed to the presence of turbulent boundary layer on 
the surface of the cylinder and the deviation at lower Re to the natural convection effects 

Ktu-niadakis [29] has studied the unsteady, forced convection heat transfer from 
an isolated cylinder up to Re =200 by numerical simulation of the Navier-Stokes and 
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energy cc] nations using spectral element method. Simulations were carried out for 
constant heat ilux and constant temperature case with the flow started impulsively. The 

variation of the length of recirculatrion zone, St, Cd with Re and the variation of 
local Nii over the surface of the cylinder were presented. 

Figueiredo et al. [30] have studied the external forced convection from an 
isothermal circular cylinder, placed at various heights above a plane surface in the range 
of Re = 1 4x 1 0^ to 40 X 1 0 \ They have presented the results for the distribution of pressure 
coefneient around the cylinder, the variation of drag and lift coefficients, the distribution 

of local Nti around the cylinder and the variation of Nu with Re for various cylinder to 
[ilatc dislances (M). It was concluded that the presence of a plane boundary introduces an 
asynnneliy on the pressure and local heat transfer distribution around the cylinder. Both 

drag and Nu decrease with decreasing cylinder to plate distance, and the influence of the 
presence of the wall on the values of drag and Nu was significant only when H/D<0.4. 

Sliin et al. [31] have studied the steady forced convection heat transfer from a 
circular cylinder in unconfmed domain to propose a unified correlation that is valid for 
heat transfer from hot cylinder to cold fluid and hot fluid to cold cylinder for Re < 40. 

They have presented the results for the variation of Nu with Re for various temperature 
loadings (Tw =900-300 K and Ta =200-300 K) for cold as well as hot cylinder. The 
variation of fluid properties significantly affected the heat transfer rate to or from the 
cylinder. As the difference between Twand Ta magnifies, wide scatter of A’m as a function 
of Re is observed due to the property variations. They have also presented a unified 
correlation for Nu as a function of Tw /Ta and verified the results for various cylinder 
diameters and various free stream velocities. 

Lange et al. [32] have studied momentum and heat transfer from a heated cylinder 
in unconfined domain in the range of Re = 0.0001 to 200 for temperature loadings (Tw /Ta 
) of 1.003-1.5. In this study, in view of very low Re, the diameters of the cylinder and 
temperature difference were chosen so as to neglect viscous dissipation and natural 
convection elTects. 'fhe results for the variation of Nu, St and Cd with Re, with the 
properties of the fluid evaluated at Tw ,Tf and Tm were presented. Drag was found to 
increase with the increase in temperature loading (Tw/ Too) with the fluid properties 
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evaluated at Iree steam temperature. The increase in drag was found to be significant at 
lower At.'. When the iluid pioperties were evaluated at cylinder wall temperature (Tw) the 
incitasc in drag was I'ound to be significant only when Re >0.1. The effect of temperature 
loading on St was found to be insignificant. In a later study by Shi et al. [33] has 
suggested that indeed there is an effect on temperature loading and property variation on 
St. 1 he transient density term that was neglected in the previous study of Lange et al [32] 
was included and the effect of temperature loading, density variation, viscosity variation 
on yVw, Ct), and St were presented for Re from 0.001 to 175. The effect of density and 
viscosity variation with temperature on Nu, Co, and St was found to be significant with 
more contribution coming from viscosity variation. The density variation showed an 
apposite el feet in steady and unsteady flows and the effect of property variation was 
signilicant only at low Re < 10. The separate characteristic temperatures for momentum 
and heal transfer were identified. 

Lange et al. [34] have extended their study to investigate the effect of the 
presence of wall on heat transfer from the probe of Hot-wire anemometer. In this study 
the range of Ae used is 0.00 1-0.1. 7 o simulate the flow conditions near wall, Couette flow 
velocity distribution is specified at the inlet. The top boundary (wall) was kept at a large 
distance so as to have no effect on the heat transfer from cylinder. On the bottom 
boundary (wall closer to cylinder), two types of boundary conditions were specified: a 
highly conducting wall and an insulating wall, to study the effect of the presence of 
conducting and non-conducting wall on the heat loss from cylinder. Wall distances (H/D) 
between 10 and 300 are considered. In the case of highly conducting wall the temperature 
gradients were high due to the proximity of highly conducting wall. Hence there was an 
increase in Nu as the cylinder approached wall. In the case of non-conducting wall, there 
was a reduction of heat transfer from the cylinder to the wall due to the suppression of 
diffusive (lux. I lence there was a decrease in Nu. 

Shi et al. [35] have studied the effect of thermal conductivity {K) of the wall and 
the shear rate on the heat loss from the cylinder for various wall distances. They have 
concluded that the presence of wall significantly affects heat transfer rate even when the 
thermal conductivily of the wall is very low. The wall effect becomes insignificant above 
1 I/I ) ■ 5. 1 lowcver, the effect ol' shear rate becomes significant. 
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Nakamura et al. [36 1 have studied the variation of Nu with flow regimes behind 
the circular cylinder experimentally. The authors have presented correlations for Nu and 
Nu at the rear .stagnation point for different flow regimes behind the cylinder; laminar 
shedding regime (Re < 150), three dimensional flow regime (300< Re < 1500) and shear 
layer transition region (3xl0^< Re < 15x10^). Nur (at rear stagnation point) followed the 
trend of base suction pressure coefficient, which depends on the flow conditions behind 
the cylinder. The exponent of Re in the correlation for mean Nusselt number for three 
dimensional How regime was lower than that for the other regimes. This is in accordance 
with the decrease in local Nu at rear stagnation point. 

Sanitijai el al. [37] have carried out experimental investigation on heat transfer 
from a circular cylinder to air and ethylene glycol/water mixtures in the range of 2xl0^< 
Re < 9x 1 0^ and 2x 1 0^< Pr < 9x10'*. The authors have considered flow regimes behind the 
cylinder and the contribution of heat transfer from front and rear parts of the cylinder 
separately a.s the flow structure is different for the two regions. They have concluded 
from their study that the three regions; 0°< 0 < 85°(growth of laminar boundary layer), 
85**< $ < 135‘*(reattachment of free shear layer), 135°< 0< 180°(the periodic vortices), 
contribute separately to the overall heat transfer from the cylinder. The heat transfer 
variation in the region of 135"< ^ < 180° also depends on the flow regime behind the 
cylinder. The authors have presented correlations for the local and overall heat transfer 
coefficients. 

Ahmed and Yovanovich [38] have obtained approximate analytical solutions for 
forced convection heat transfer from circular cylinders, flat plates and spheres in cross 
flow, 'fhey have proposed a method to predict heat transfer that is obtained by 
lincari/ation o(' energy equation and then transforming the linearized energy equation into 
transient heat conduction equation, which has analytical solutions. The analytical 
solutions were presented in closed form for wide range of Re and Pr and results were in 
good agreement with the previous experimental results and analytical solutions. 

Similarly, Khan et al [39] have investigated the effect of blockage ratio on fluid 
flow and heat transfer from a circular cylinder confined in a plane channel for isothermal 
and isoflux boundary conditions using analytical approach. Correlations were presented 
for drag coefficient as a function of blockage ratio and Re, and for heat transfer 
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coefficient as a function of Re, Pr and blockage ratio. Separate correlations were 
presented for isothermal and isoflux boundary conditions. Nu and drag were found to 
increase with increase in blockage ratio and Nu for isoflux boundary condition was found 
to be higher than Nu for isothermal boundary condition for all blockage ratios. 

In light of the above review, it may be concluded that there exist few studies in 
the literature for the momentum transfer from an asymmetrically confined circular 
cylinder [24] and heat transfer from a symmetrically confined cylinder [39] in a plane 
channel. There have been no studies in the literature regarding the heat transfer from an 
asymmetrically confined cylinder in a plane channel. Therefore, to fill the existing gap, 
we have studied in the present work the momentum and heat transfer from an 
asymmetrically confined circular cylinder in a plane channel using FLUENT 6.1. 
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CHAPTERS 

PROBLEM STATEMENT AND FORMULATION 


3.1 Two-dimensional iso-thermal flow past a circular cylinder in unbounded 
domain: 

Consider the unsteady 2-D isothermal flow of an incompressible Newtonian fluid 
past a circular cylinder in unbounded domain, as shown schematically in Figure 3.1. The 
geometry for 2-D flow in unbounded domain consists of a circular cylinder of radius r, 
placed in a circular outer boundary of radius Rco, where R«/r = 140. The radius of outer 
circular boundary (i?®) is taken to be large enough to ensure that the effect of boundaries 
on the inner cylinder wall is negligible. The left half of the outer circular boundary is 
taken as velocity inlet and the right half as the outflow. 

The objectives of this 2-D study are: (a) to study the transition of flow past a 
cylinder from steady to unsteady flow, characterized by Von karman vortex street, with 
the aid of various upwinding schemes available in FLUENT, and (b) to determine the 
critical Re at which the flow becomes unsteady using each of the upwinding schemes. 

3.2 Three-dimensional iso-thermal flow past a circular cylinder in unbounded 
domain: 

The geometry taken to simulate 3-D isothermal flow past a circular cylinder in 
unbounded domain, as shown schematically in Figure 3.2, is a rectangular brick of 
dimensions 28c/, \6d and 125d in x, y and z-directions, respectively, where d is the 
diameter of the cylinder. The cylinder is located at lOc/and 18c/ from the inlet and outlet, 
respectively. The fluid enters at the plane ABCD with a uniform velocity U® and leaves at 
the plane EFGIL 

I'he main objective of this 3-D study is to examine the transition of flow past a 
cylinder from 2-D unsteady to 3-D unsteady flow, characterized by the dislocation of 
alternating vortices i.e. the bending of vortices in z-direction. The study of heat transfer 
effects is excluded in 2-D and 3-D study. 
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3.3 Non-iJvOtlicrin^I flow psist & circulsr cyiindcr conlHincd in sl plnnc chnnnciz 

In the cdse ol non-isothermal flow past a circular cylinder confined in a plane 
channel, as shown schematically in Figure 3.3, two plane walls separated by a distance H 
confine a circular cylinder of diameter d whose position is defined by gap ratio, y, defined 
as AliH/2-D/2) where A is the minimum distance from the surface of the cylinder to the 
nearest wall. 1 he value of y is equal to I when the cylinder is placed symmetrically 
between the plane walls and is equal to zero when the cylinder touches one of the walls. 

1 he cylinder i.s located at ISJ and 40<i from the inlet and outlet, respectively. The length 
ol the cylinder in ^-direction is assumed to be long enough to have insignifieant end 
elfects, implying that there is no flow in z-direction and no flow variables depend upon z- 
coordinatc. The lluid (Air, Pr = 0.7441) at a temperature of T* enters at the velocity inlet 
boundary AC with a parabolic velocity profile U = Urmjx{C(ylHI2)^) and leaves at the 
outflow boundary ED. The confining walls are at free stream temperature (T*) whereas 
the cylinder is at a temperature of Tw The effect of variation in density and viscosity is 
neglected owing to the small temperature difference (Tw - T«, = 2° K) used in the study. 

In this study, H and yare varied to study the wall effects on momentum and heat 
transfer from the cylinder at a constant wall temperature. With respect to momentum 
transfer, the transition of flow past a cylinder in plane channel from 2-D steady to 2-D 
unsteady flow is studied. The critical Re at which the flow becomes unsteady resulting in 
periodic vortex shedding is also determined for various blockage ratios d/H) and y. 
In addition to this, the effect of /? and yon drag coefficient and St is studied. With respect 
to heat transfer from the cylinder, the effects ofy^and yonNuw are determined. 

3.4 Governing Equations 

A Cartc.sian coordinate system is chosen. The velocity components in x, y and z 
directions are V^, Vy and V,. respectively. The velocity components in x-, y-, z-directions 
for 2-1 ) flow arc 

Vx =" Vx (x, y), Vy = Vy (x, y) and V^ = 0 and the velocity components for 3-D flow are 
given as 

Vx Vx (x, y, z), Vy = Vy (X, y, z) and V, = Vz (x, y, z). 
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Under these conditions the continuity equation for an incompressible fluid in 3-D is given 
as 


dv^. d^y dF. 

+ --r--’ 4- = 0 


dx 




dz 


(3.1) 


and, the corresponding x-, y- and z-components of the Navier-Stokes equations (for 
constant p and jJ) are written as: 
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Substitution of V^, = 0 in the above equations results in 2-D governing equations, 
for uon-isothermal How over the cylinder in plane channel, the 2-D energy equation is 
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3.5 Boundary conditions: 

The following are the boundary conditions for three cases, 2-D isothermal flow in 
unbounded domain, 3-D isothermal How in unbounded domain and 2-D non-isothermai 
flow in a plane channel. 


3.5. 1 Two-dinicnsiomil (2-D) isothermal flow over a cylinder in unconfmed domain: 
Inlet hotmdary: The fluid is assumed to enter the left half of the outer circular boundary 
at a uniform velocity Uoo 

Vk = U«; Vy = 0 (3-6) 

Outflow boundary: the right half of the outer circular boundary is taken as outflow 
boundary. 
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In tills case zero diffusion flux for all variables is assumed. This implies that the 
conditions of the outflow plane are extrapolated from within the domain and have no 
impact on the upstream flow. 1 his is in consistent with fully developed flow assumption. 
Surface of (he cylinder. At the surface of the cylinder, no slip boundary condition is 
applied. 

Vx = 0 ; Fy = 0 (2.7) 


3.5.2 I hrce~dimensional (3-D) isothermal flow over a cylinder in unconfined domain: 
Inlet boundary: 'fhe fluid enters the plane ABCD at a uniform velocity Uoo 
V, lU Vy -Oand V, = 0 (3.8) 

Free .slip boundaries: At the planes confining the cylinder in y-direction, ABEF and 
(IKllh free slip boundary condition is specified. This boundary condition ensures that 
the change of flow variables at the boundary in y-direction is negligible. 


9.1tj 
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(3.9) 


Symmetric boundaries: At the planes confining the cylinder in z-direction, BFGC and 
AEHD., symmetric boundary condition is specified. FLUENT assumes a zero flux of all 
quantities across a symmetry boundary. There is no convective flux across a symmetry 
plane, i.e. the normal velocity component at the symmetry plane is zero. There is no 
diffusion flux across a symmetry plane, i.e. the normal gradients of all flow variables are 
zero at the symmetry plane. The symmetry boundary condition therefore ensures zero 
normal velocity and zero normal gradients of all variables at the symmetry plane. 

Outflow boundary: At the plane EFGH, outflow boundary condition is specified. 

Surface of (he cylinder: At the surface of the cylinder, no slip boundary condition is 
applied. 

0; f;, = 0 and V, = 0 (3-10) 


3.5.3 'I'wo-dimcnsional (2-D) non-isothermal flow in a plane channel. 

Met boundary at AC: 'fhe fluid is assumed to enter the inlet with a parabolic velocity 
prof! le { f ( 1 -(y////2) ^) with a free stream temperature of Too. 

Uv ^ I fun ( 1 -(y////2) ~); Fv = 0 and T = T« (3-11) 
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Outflow boundary at BD: Zero diffusion flux for all variables is assumed at the outlet 
boundary. 

Walls AB and C.£). No slip condition with impervious wall is applied at confining walls 
AB and 

CD and these walls are at free stream temperature of T*. 

-O; =0andT = Too (3.12) 

Surface of the cylinder: At the surface of the cylinder, no slip boundary condition is 
applied and tlic cylinder is at a temperature of Tw 

= 0 andT^Tw (3.13) 


3.6 C’alcuhition of drag and lift coefficients 

I'hc total drag force exerted by the fluid on the cylinder is the sum of pressure and 
friction drags, flic total drag coefficient (dimensionless) in case of unbounded domain is 
defined as 


/ipuid 


(3.14) , 


where Fpis the total drag force per unit length of the cylinder. 

The lift coefficient (dimensionless) in case of unbounded domain is defined as 
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(3.15) 


where Fo is the lift force per unit length of the cylinder. 

Tlie total drag coefficient in case of a cylinder confined in a plane channel is defined as 


(3.16) 

where is tlic average velocity based on the diameter of the cylinder given by 



where d/ and d 2 arc y-coordinates of the cylinder and 
below 


(3.17) 
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(3.18) 


rf, Ml- ,») -!)+«' 


(3.19) 


The lift coefficient (dimensionless) in case of confined domain is defined as 


Q = 


YifFild 


(3.20) 


3,7 Calculation of Strouhal number 

The non-dimensional frequency of oscillation of cylinder is defined as Strouhal number. 
In case of unbounded domain, Strouhal number is given as 


til 


/ d 


(3.21) 


Strouhal number in the case of confined domain is given as 


57 


/ 

vjd 


(3.22) 


where /is the frequency of oscillation of the lift coefficient of the cylinder or the vortex 
shedding frequency. 
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Figure 3.1: Schematic for 2-D flow past a cylinder in unbounded domain (RJr = 140) 
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Figure 3.3: Schematic for flow past an asymmetrically confined cylinder in a plane channel 
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Chapter 4 

Numerical Methodology 


In the present study, a commercial flow modeling software, FLUENT 6.1 has 
been used to cany out numerical simulations. FLUENT employs a finite volume based 
numerical method to solve the governing flow and energy equations. 

4.1 Finite volninc method 

I he finite volume based method consists of (a) division of the domain into 
di.screte control volumes using a computational grid, (b) integration of the governing 
cc|uations on the individual control volumes to construct algebraic equations for the 
discrete de|icndcnt variables such as velocities, pressure, temperature, and (c) 
linearization of the discretized equations and solution of the resultant linear equation 
system to yield updated values of the dependent variables. 

4.2 Grid generation 

The first step in the solution of either flow or heat transfer problem using 
FLUENT is to build the appropriate geometry and divide the geometry into discrete 
control volumes (grid or cells). In this study, a commercial grid generation software, 
GAMBIT (in built with FLUENT), is used for the creation of geometry, division of the 
domain into discrete control volumes, and for the specification of the appropriate 
boundary conditions on the edges (or faces). 

4.2. 1 Two-dhmn.sional grid generation 

The grid generation procedure for 2-D geometries is summarized as follows: 

• The edges arc created from the vertices followed by faces from edges. 

• 'fhe edges arc graded using different grading schemes available in GAMBIT, namely: 
succes.sive ratio, first last ratio, and last first ratio. 

• fhe faces (the computational domain) are meshed using different paving schemes 
avaihible in GAMBIT, namely: Quadrilateral-Pave and Triangular-Pave scheme. 

• Appropriate boundary conditions are specified on the boundaries (edges). 
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I [it computdtional domain in the case of rectangular geometry is divided into 
discitlt contiol volumes using either ot the schemes Quadrilateral-Pave or Triangular- 
Pavc scheme. I he combination ot these two schemes, Quad-Tri-Pave scheme is used for 
2-1) circular geometry. The accuracy of any numerical solution depends on the fineness 
of the grid used for the computation. A non-uniform mesh with finer mesh around the 
cylinder where the velocity gradients are expected to be high and coarse mesh away from 
the cylinder is created using successive ratio grading scheme on the edges. The grid 
sensitivity analysis is carried out by increasing the number of uniform intervals on the 
surlace o( the cylinder until the drag coefficient and Nusselt number became insensitive 
to further mesh refinement. 

4.2.2 Threc-dimcmkmal grid generation 

The grid generation procedure for 3-D geometries is summarized as follows: 

• The volumes are created using readily available geometries such as a cylinder and a 
rectangular brick. Boolean operations are performed on the volumes to define 
appropriate computational domain. 

• The volume is divided into smaller volumes so as to facilitate finer mesh around the 
cylinder and coarser mesh away from the cylinder. 

• The edges are graded using different grading schemes available in GAMBIT. 

• The faces are meshed using Quadrilateral-Pave schemes followed by meshing the 
volumes using Cooper scheme which creates identical mesh in the third (z) direction. 

• Appropriate boundary conditions are specified on the boundaries (faces). 

4.3 Solution procedure 

After specifying the proper boundary conditions on the edges or faces, the mesh 
file is exported to FLUENT for simulation purposes. The solution method involves a few 
steps as outlined below. 

4. 3. 1 Selection of the solver 

FLUENT allows user to choose either of the two numerical solution methods, 
namely: (a) segregated solver, and (b) coupled solver, to linearize and solve the 
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discietized governing equations. Both of the solution methods use a finite volume 
disci etization method. However, the approach to linearize and solve the discretized 
equations is different for both the schemes. In general both the solvers are applied in 
solving wide variety of fluid flow and heat transfer problems. The segregated solver is 
specifically used in solving incompressible flows and mildly compressible flows, 
whereas the coupled solver is specially developed for solving high-speed compressible 
flows. However, in this study, keeping in view of the incompressible nature of the flow, 
the segregated solver is used. In this method, the governing equations, namely: the 
continuity, momentum and energy equations (heat transfer) are solved sequentially i.e. 
one after the other (one at a time). 

4.3,2 Selection of the ditscretization scheme 

4.3.2. 1 Pre.ssure interpolation scheme 

FLUliN f uses a collocated grid in which both vectors (velocity and fluxes) and 
scalar quantities (pressure and temperature) are stored at cell centers. However, in the 
discretized momentum equation, face values of pressure are required. Therefore, an 
interpolation scheme for pressure is required to compute the face values of from the 
adjoining cell values. I'here are five interpolation schemes available in FLUENT for 
pressure computation: (i) standard, (ii) linear, (iii) second order, (iv) body-Force 
weighted and (v) PRESTO. In the present study standard pressure interpolation scheme 
has been used. 

4. 3. 2. 2 Pressure-Velocity coupling 

In segregated solution method, the continuity equation is used as an equation for 
pressure, even though pressure does not appear explicitly in the discretized continuity 
equation. Therefore some numerical scheme is required to account for pressure velocity 
coupling. There are three algorithms available in FLUENT for pressure-velocity 
coupling: (i) SIMPLE, (ii) SIMPLEC, (iii) PISO. The pressure-velocity coupling scheme 
pressure-implicit with splitting of operators (PISO) is used for unsteady flows as 
suggested by FLUENT user manual. 
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43.2.3 Momentum Upwinding 

As mentioned earlier, FLUENT uses a collocated grid in which both fluxes and 
scalar quantities are stored at cell centers. However, for convective terms in the 
discretised momentum and energy equations, the face values of the dependent variables 
are required which are interpolated from the cell center values of the adjoining cells. 
Upwinding means that the face values are derived from cells in the up wind direction 
relative to the direction of the flow. There are four upwinding schemes available in 
1‘LUbNI: (i) first order upwind, (ii) second order upwind, (iii) power-law, and (iv) 
QUICK. In this study second order upwinding scheme is used for both the momentum 
and energy equations. In this method the cell face values is approximated by expanding 
through a 'faylor series about the cell center value. The comparison of three-upwinding 
schemes, first order, second order and QUICK is also studied with respect to the 
transition of flow past a circular cylinder from steady to unsteady flow. 

433 Convergence criteria 

In the present study, a convergence criterion of lx 10"°^ has been applied to the 
residuals of the continuity, the momentum and energy equation. In the case of steady 
flows, the iterations are stopped when the solution reaches steady state, whereas in the 
case of unsteady flows, the iterations are stopped when the oscillations of the lift 
coefficient reach a periodic steady state. 
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CHAPTER 5 

MODEL VALIDATION 

In the present study, the momentum and heat transfer from a circular cylinder has 
been studied using the commercial flow modeling software, FLUENT 6.1. Before solving 
the problem of interest, the code has been validated with the results reported in the 
literature. A reliable validation helps in determining the accuracy of the code used and 
gives the necessary idea about the problem to be solved. In this chapter, the comparison 
of our results with the published data for the steady and unsteady momentum and heat 
transfer from a circular cylinder in confined as well as in unconfined domain is presented. 

5.1 Momentum and heat transfer from a cylinder in unconflned domain in the 
steady flow range {Re < 50) 

For the validation of our results with the published data for Re< 50, a circular 
geometry with a circular outer boundary located at 1 12 radii way from the surface of the 
cylinder is chosen. A grid size of 90x60 previously used by Chakxaborty et.al [55] is 
chosen for the computation in this work. The values of Cd obtained in this study are 
compared with the published results and are found to be in good agreement, as shown in 
Table 5.1 and Figure 5.1. The maximum deviation from the literature values is 8.2% at 
Re = 40, which is due to the smaller outer boundary ( xjx = 29) used by Sucker and 
Brauer [42] 

The validation of heat transfer results for Nuw { Pr = 0.7441) for the flow of air 
past a circular cylinder in unconfined domain is carried out using a circular geometry, 
with the outer circular boundary fixed at 1000 radii away from the circular cylinder, 

keeping in view of the low Re. The values of Nuw for the range of Re between 0.01 and 
40 are compared with the published results and are found to be in good agreement as 
shown in Table 5.2 and Figure 5.2. A maximum deviation of 4.1% from the literature 
value is observed at Re = 40, which is within the acceptable accuracy. 
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5.2 Momentum transfer from a cylinder in unconfined domain in the unsteady flow 
range 

The validation of Cj, (time average). Cl and St for the range of Re from 50 to 
1 000 is carried out using a circular geometry, with the outer circular boundary located at 
140 radii away from the cylinder. The grid sensitivity analysis is carried out on this 
geometry, as shown in Table 5.3. It is observed that as the grid size is refined from 
1 10x120 to 120x144, the change in the value of is less than 0.24%. Hence, the 
values obtained using the grid size of 120x144 may be treated as grid independent. 
Therefore, the grid size of 120x144 was used as the optimum grid for further 
computations in this study. In Table 5.4 and Figure 5.3, values computed in this study 
are compared with those available in the literature. The maximum deviation in C^, from 
the literature values for Re below 500 is 5.3%. At Re =1000, the deviation in value 
from that of Braza et al. [17] is 23%. However, the drag value at Re =1000 is in good 
agreement with values obtained by Franke et al. [22] and Henderson [40], with a 
deviation of less than 2.4%. The reason for the large deviation (23%) of value 
computed by Braza et al. [17] from that of the values reported in the literature is not 
known. 

In Table 5.5, Cl and St values are compared with the values reported in the 
literature. As observed, Cl is in good agreement with the values obtained by Franke et al 
[22]. The comparison of Cl with those obtained by Braza et al. [17] could not be made as 
the authors have reported r.m.s values of Q. It is observed from Table 5.5 that there is a 
good agreement between the values of St calculated in the present study and the literature 
at all Re, except at Re = 1000, at which the calculated value is 7.7% higher than the value 
obtained by Braza et al. [17]. The comparison of St values with those obtained by Braza 
et al. [17] and Franke et al. [22] is also shown in Figure 5.4. 

The comparison of , Cl and St for Re = 50, 1 00, 200, and 1 000 with various 
simulation and experimental investigations reported in literature is presented in Table 5.6. 
Cd value obtained in this study for Re = 50 is in excellent agreement with the simulation 
results of Saiki and Biringen [43] and that of the experimental result of Clift et al. [44], 
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with a maximum deviation of 1 .2%. However, at Re = 50 there is a deviation of 23% in 
value that of Gresho et al. [21]. It is to be noted from Table 5.6 that obtained by 
Gresho et.al. [21] for =100 and 200 also show a large deviation from the rest of the 
values reported in the literature, whereas obtained in our study is in good agreement 
(deviation < 2.05%) with the experimental value of Clift et al. [44] as well as with the 
simulation results of Kim et al. [47], and Wanderley and Levi [48]. reported in the 
literature at Re =200 is in the range of 1.24 to 1.76, as shown in the Table 5.6 and 
obtained in our study is 1 .298 1 , which is within the range. 

From Table 5.6, St value for Re = 50 is observed to be in good agreement with the 
experimental values of Roshko [45], and Berger and Wille [46]. In the present study, the 
maximum deviation of St value from experimental value at Re = 50 is less than 1.6%. The 
value of St at =100 is in good agreement with the experimental values of Roshko [45], 
and Berger and Wille [46], as well as with the simulation results presented in Table 5.6. 
The maximum deviation at Re =100 is 7.5% from the simulation result of Rengel and 
Sphaier [49]. The value of St at Re = 200 is in good agreement with the experimental as 
well as simulation results, with a maximum deviation of 9% from the value obtained by 
Gresho et al. [21]. Finally, C^and St values at Re =1000 are found to be in good 
agreement with the literature values, with a maximum deviation of 2.4% and 5.3%, 
respectively. 

5.3 Momentum and heat transfer from a cylinder confined in a plane channel 

In the present thesis, the main objective is to study the momentum and heat transfer 

characteristics (C^, St, and Nuw) from an asymmetrically confined circular cylinder in a 
plane channel. Before solving the problem of our interest, the results for momentum and 
heat transfer for an asymmetrically confined circular cylinder in plane channel are validated 
with the published data of Zovatto and Pedrizzetti [24] for P= 0.2 and y= 0.375 and also 1 
(symmetrical case). The simulation of flow past a symmetrically confined as well as 
asymmetrically confined cylinder is carried out using the rectangular geometry with a 
downstream length of 40^/ and upstream length of \5d, and with a parabolic velocity 
specified at the inlet, as shown in Figure 3.3. Referring Table 5.9 for 0.2 and y = 0.375 
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the drag values are found to be in good agreement with the results obtained by Zovatto and 
Pedrizzetti [24]. The maximum deviation in Cj, is 5.8% at Re = 40, which is within the 
acceptable accuracy. 

Since the published data available in the literature for an asymmetrically confined 
cylinder is limited, the model was also validated with the data of Khan et al [41] for a 
symmetrically confined cylinder in a plane channel. The comparison of results for and 

Nuv with the published data of Khan et al. [41] for /?= 0.1 and 0.2 is shown in Tables 5.7 
and 5.8, respectively. The comparison is also described in Figures 5.5 and 5.6. The 
deviation in C^, for >5= 0.1 obtained in this study from those obtained by Khan et al. [41] 

is significantly high with a maximum deviation of 40% at = 50, whereas Nu^,. values 
are in close agreement \vith a maximum deviation of 4.7% at Re = 20. The maximum 

deviation in for y6= 0.2 is 24% at Re = 69, whereas the maximum deviation in Nuy,. is 
9.5% at Re =1000. This large deviation in is attributed to the difference in the velocity 
profile specified at the inlet of the channel and different velocity profile neair the cylinder. 
In the study of Khan et al. [41], a uniform velocity is specified at the inlet and near the 
cylinder, whereas in our study a parabolic velocity is specified at the inlet and the cylinder 
is also assumed to experience the parabolic velocity profile. The reference velocity used for 
the drag calculation in our study is the velocity averaged over the surface of the cylinder, 
while Khan et al. [41] have used the free stream velocity as the reference value. The above 
mentioned differences in the two studies have caused deviation in . The drag values for 
yff = 0.2 are also found to be in good agreement with those obtained in the study by Zovatto 
and Pedrizzetti [24], in which case the parabolic velocity is specified at the inlet and the 
average velocity is used as the reference for the drag calculation. 

From these comparisons, it may be concluded that and 77 m obtained using 
this code are accurate within 5-10% of the values reported in the literature. The close 
agreement between the results obtained in this study and the literature gave us the 
confidence in the code used. Therefore, we have used the same code for studying the 
momentum and heat transfer characteristics from an asymmetrically confined circular 
cylinder in a plane channel for various blockage ratios (Ji) and gap ratios (y). 
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S.No 

Re 

Cd 

This 

work 

Cd 

Chakraboity et 

al. 

(2004)* 

Cd 

Kawaguti and 

Jain 

(1966)’* 

Cd 

Sucker and 

Brauer 

(1975)^ 

Cd 

D’ Alessio 

and Dennis 

(1994)^ 

1 

1 

10.374 

10.445 

- 

- 

“ 

2 

10 

2.757 

2.803 

2.870 

- 

- 

H 

20 

2.001 

2.035 

2.090 

2.178 

1.941 

■■ 

30 

1.702 

1.702 

1.760 

“ 

- 

5 

40 

1.499 

1.517 

1.529 

1.633 

1.443 

6 

50 

1.392 

1.392 

1.440 


- 


Table 5.1: Comparison of Cd values with the published data for the steady flow 
range ( rojr = 112, ** VaJr = 112, ^ taJr = 29, r^Jr = 55) 


S.No 

Re 

Nuw 

This work 

Nuw 

Dennis and Hudson 
(1968) 

Nuw 

Collis and Williams 
(1959) 

Nuw 

Lange et al. 
(1998) 

1 

0.01 

0.308 

0.302 

0.296 

0.29 

2 

0.1 

0.453 

0.452 

0.439 

0.44 

3 

1 

0.819 

0.812 

- 

0.80 

4 



1.897 

" 

1.9 

5 


mg^mi 

2.557 

- 

2.5 

6 



3.480 

- 

3.4 


Table 5.2: Comparison of Nuw (Pr = 0.7441) values with the published data for the 
steady flow range 


S.No 

Re 

Grid size (90x120) 

Grid size (11 Ox 120) 
Cd 

Grid size (120x144) 

1 

20 

1.956 

1.949 

1.948 

2 

40 

1.513 

1.513 

1.512 

3 

50 

1.390 

1.393 

1.397 

4 

100 

1.302 

1.302 

1.305 


Table 5.3: Grid independence test using QUICK scheme for momentum upwinding 
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S.No 

Re 

Co 

This work 

Co 

Braza et al. 
(1986) 

Co 

Franke et al. 
(1990) 

Henderson 

(1995) 

Max % deviation 
with literature 
value 

1 

50 

. 1.393 

1.400 

1.39 

1.412 

1.3% 

2 

100 

1.302 

1.253 

- 

1.349 

3.9% 

3 

150 

1.289 

- 

- 

1.333 

3.3% 

4 

180 

1.293 

- 

- 

1.336 

3.3% 

5 

200 

1.298 

1.321 

1.31 

1.341 

3.2% 

6 

250 

1.314 

1.359 


1.357 

3.2% 

7 

300 

1.329 

1.405 

1.32 

1.376 

5.3% 

OO 


1.381 



1.448 

4.5% 

9 

1000 

1.473 

1.198 

1.47 

1.509 

23% 


Table 5.4: Comparison of values with the published data for the unsteady flow 
range 


S.No 

Re 

Cl 

This work 

Cl 

Franke et al. 
(1990) 

St 

This work 

St 

Braza et al. 
(1986) 

St 

Fraitke et al. 
(1990) 

1 

50 

±0.010 

±0.0 

0.1212 

0.120 

0.116 

2 

100 

±0.269 

- 

0.1600 

0.165 

- 

3 

150 

±0.455 

- 

0.1790 

0.181 

- 

4 

180 

±0.537 

- 

0.1870 

0.185 

- 

5 

200 

±0.599 

±0.65 

0.1912 

0.198 

0.194 

6 

250 

±0.701 

- 

0.1990 

0.204 

- 

7 

300 

±0.814 

±0.84 

0.2060 

0.209 

0.205 

8 

500 


1 

0.2222 

- 

- 

9 

1000 

±1.307 

±1.36 

0.2370 

0.22 

0.236 


Table 5.5: Comparison of Cl, 5/ values with the published data for the unsteady 
flow range 
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Re 

Cz, 

Cl 

St 

Tseng and Ferziger (2003) 

100 

1.42 

0.29 

0.164 

Kim et al. (2001) 

100 

1.33 

- 

0.165 

Saiki and Biringen (1996) 

50 

1.38 


0.139 


100 

1.26 

- 

0.171 


200 

1.18 

- 

0.192 

Gresho et al. (1984) 

50 

1.81 

- 

0.14 


100 

1.76 

- 

0.18 


200 

1.76 

- 

0.21 

Clift etal.(1978) 

50 

1.41 

- 

- 

(Experimental) 

100 

1.33 

- 

- 


200 

1.24 

- 

- 

Roshko (1953) 

50 

- 


0.122 

(Experimental) 

100 

- 

- 

0.167 

Berger and Wille (1972) 

50 

- 

- 

0.12-0.13 

(Experimental) 

100 

- 

- 

0.16-0.17 


200 

- 

- 

0.18-0.19 

Meneghnini (1993) 

100 

1.520 

0.353 

0.162 


200 

1.395 

0.570 

0.195 

Herfjord (1995) 

100 

1.36 

0.34 

0.168 


200 

1.35 

0.70 

0.196 


1000 

1.47 

1.45 

0.234 

Rengel and Sphaier (1999) 

100 

1.36 

0.32 

0.173 


200 

1.35 

0.67 

0.203 


1000 

1.50 

1.70 

0.225 

Wanderley and Levi (2002) 

100 

1.33 

0.195 

0.163 


200 

1.32 

0.534 

0.192 


1000 

1.51 

1.370 

0.235 

Present study 

50 

1.393 

0.010 

0.121 


100 

1.302 

0.269 

0.160 


200 

1.298 

0.599 

0.191 


1000 

1.473 

1.307 

0.237 


Table 5.6: Comparison of C,, , Cl, St values with the experimental and simulation 
results 
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S.No 

Re 

C. 

This work 

Co 

Khan et al. (2004) 

Niiw 

This work 

Nil. 

Khan et al. (2004) 

1 

20 

2.157 

2.526 

2.526 

2.411 

2 

30 

1.780 

2.266 

2.969 

2.953 

3 

40 

1.567 

2.112 

3337 

3.410 

4 

50 

1.430 

2.009 

3.660 

3.813 

5 

100 

1.296 

1.754 

5.219 

5.392 

6 

200 

1.290 

1.577 

7.596 

1.626 

7 



1.499 

9.493 

9.340 

8 

400 

1.338 

1.453 

11.01 

10.78 

9 

500 

1.340 

1.422 

12.40 

12.05 

10 

1000 

1.40 

1.344 

17.30 

17.05 


Table 5.7: Comparison of C^, and Nu^ (Pr = 0.7441) values for 0.1 


S.No 

Re 

Co 

This 

work 

Co 

Zovatto and 
Pedrizzetti (2001) 

Co 

Khan et al. 
(2004) 

Nuw 

This work 

Niu- 

Khan et al. 
(2004) 

1 

20 

2.671 

- 

2.603 

2.667 

2.450 

2 

30 

2.142 

- 

2.337 

3.119 

3.000 

3 

40 

1.860 

- 

2.180 

3.497 

3.464 

4 

50 

1.681 

- 

2.074 

3.830 

3.873 

5 

60 

1.556 

- 

1.996 

4.135 

4.243 

6 

69 

1.471 

1.6 

1.941 

4.4112 

4.550 

7 

200 

1.253 

1.3 

1.632 

7.592 

7.747 

8 

300 

1.268 

1.2 

1.553 

9.326 : 

9.489 

9 

400 

1.274 

- 

1.505 

10.91 

10.95 

10 

500 

1.324 

- 

1.473 

12.13 

12.21 

11 

1000 

1.368 

- 

1.395 

18.98 

17.32 


Table 5.8: Comparison of and Nuv (Pr = 0.7441) values for y^= 0.2 


35 

















S.No 

Re 

Co 

This work 

Zovatto and Pedri 2 zetti 
(2001) 

Nuw 

This work 

1 


3.150 

3.0 

2.706 . 

2 

40 

2.222 

2.1 

3.609 

3 


1.852 

1.8 

4.269 

4 

108 

1.481 

1.5 

5.539 

5 

200 

1.333 

1.3 

7.615 


Table 5.9: Comparison of and Nu^ (iV = 0.7441) values for and ^=0.375 
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CHAPTER 6 

RESULTS AND DISCUSSION 

This chapter presents the simulation results for the flow of an incompressible, Newtonian 
fluid past a circular cylinder in unconfined as well as confined domains. The main objective of 
this study, as mentioned earlier, is to study the transition of flow fi-om 2-D steady to 2-D 
unsteady in unconfmed and confined domains. In the case of flow in unconfined domain, the 
validation of simulation results for Cj, , St, and critical Re with the published data is presented. 
In addition, the comparison of and St and critical Re obtained using various upwinding 
schemes available in FLUENT is presented. Since the data available in the literature for the 
transition of the flow past a circular cylinder confined in a plane channel is limited [23,24], the 
transition is characterized as the function of asymmetry ( 7 , gap ratio) over a wide range of 
blockage ratios (/^. In this work, in particular the critical Re at which transition in the flow takes 
place is determined as a function of /and /?. 

6.1 The transition of flow past a circular cylinder from steady to unsteady in unbounded 
domain 

The flow past a circular cylinder may be divided into different regimes depending on Re. 
For Re ~ I, the flow is fore and oft symmetric, divides at front stagnation point and reunites 
smoothly behind the cylinder. As Re increases above 5, the flow separation takes places on the 
downstream side with formation of two symmetric vortices, which are stable and remain attached 
to the cylinder. Around Re = 47, flow transition takes place from steady to unsteady, in which 
case the oscillations of vortices give rise to the well-known Von karman vortex street, with the 
non-dimensional number St characterizing the frequency of oscillation of the vortices. In this 
section, this flow transition is numerically studied by the various upwinding schemes available in 
FLUENT. The critical Re at which the transition takes place is determined for each of the 
upwinding schemes by carrying out unsteady numerical simulations. 

The grid (120x144) used for the unsteady simulations is shown in Figure 6.1, which was 
found to be optimum by carrying out grid sensitivity analysis. The results for the grid 
independence study are shown in Table 5.3. In the simulation, the radius (i?«) of outer circular 
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boundary is taken to be large (KJr = 140) enough to ensure that the effect of boundaries on the 
flow field near the cylinder wall is negligible. The resolution of boundary layer near the cylinder 
wall significantly influences the numerical results. Therefore, to resolve the boundary layer 
accurately, a constant expansion factor is used while meshing the geometry that produces a non- 
uniform grid with a finer mesh near the cylinder and coarser mesh away from the cylinder. In 
particular, the distance of the first grid point fi'om the cylinder wall is kept at 0.007% of the 
diameter of the cylinder. The time step (At) in the unsteady simulation is kept below T/lOO, 
where T is the time period of the oscillations, i.e. more than 100 time steps per period are used. 

6. 1. 1 Vortex shedding initiation mechanism: 

The time marching calculations were started with the step increase in the velocity 
from zero to free stream value of Ua.. For Re > 40, one experimentally observes Von karman 
vortex street which is the characteristic of the unsteady flow. In the case of numerical 
simulations the destabilizing factors like, non-uniform inlet conditions, irregularity of the 
boundary conditions, vibrations and surface roughness that are usually present in the 
experimental conditions are absent as the boundary conditions and initial condition are 
symmetric, which in general should lead to the steady symmetric solution of the Navier-Stokes 
equations. However, in numerical simulations, the truncation errors, round-off errors and the 
errors due to numerical approximations act as destabilizing factors, which lead to vortex 
shedding. But the time period required for numerical errors to accumulate and eventually lead to 
vortex shedding turns out to be significantly large. Therefore, to reduce this transition period 
before the periodic vortex shedding could be observed fi'om the cylinder, Braza et al. [17] have 
introduced a perturbation method to generate vortex shedding fiom the cylinder, which is based 
on the artificial rotation of the circular cylinder in clockwise direction followed by rotation in 
anticlockwise direction for a short duration and then keeping the cylinder stationary for the 
remaining part of the calculations, as depicted in Figure 6.2. 

In our study, we have also used the method of Braza ef al. [17] to generate vortex 
shedding from the cylinder. It has been shown by Martinez (1979) that any random or imposed 
perturbations die out for Re < 40, whereas they result in vortex shedding for Re > 50. The final 
periodic state solution should be the same irrespective of the magnitude of the disturbance, if the 
vortex shedding is intrinsic property of the flow at Re > 50. To ascertain this, the effect of the 
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magnitude of the disturbance has been tested by varying the speed of the rotation of the cylinder 
for Re =100, as shown in Table 6.1. The cylinder is rotated for about hundred time steps in 
clockwise direction follovved by rotation in anticlockwise direction using two different rotational 
speeds, 2 rad/sec and 6 rad/sec. It is evident from the results presented in Table 6. 1 that the effect 
of the magnitude of the disturbance on global pareimeters like C ,^ , Cl and St is very 
insignificant. The final periodic state reached by the oscillations of and Cl for the two 
disturbances may be clearly observed from Figure 6.3. They are essentially identical for two 
disturbances, one of which is three times the other. From these observations, it may be concluded 
that vortex shedding is the intrinsic property of the flow and the final periodic state solution is 
not influenced by the magnitude of the disturbance introduced into the flow. 

6.1.2 Comparison of upwinding schemes: 

The comparison of various upwinding schemes available in FLUENT for momentum 
upwinding, namely; first order upwinding, second order upwinding and QUICK schemes is 
presented with respect to the transition of the flow from steady to unsteady flow. The flow is 
perturbed slightly and the time dependent motion following perturbation is examined. The 
critical Re at which the disturbance given through the rotation of the cylinder amplifies leading 
to periodic vortex shedding is determined for each of the upwinding schemes. The results are 
shown in Table 6.2. The critical Re determined using second order and Quick upwinding 
schemes is 50, which is very close to the value (46) obtained by Jackson [18] using finite 
element method. The history of lift coefficient for Re = 48 and 50 is shown in Figure 6.4 for the 
second order upwinding scheme. The disturbance introduced into the flow through the rotation of 
the cylinder decayed at Re = 48, whereas at Re = 50 it is amplified, leading to periodic vortex 
shedding from the circular cylinder. The unsteady nature of the flow at Re = 50 is depicted in 
Figure 6.5, in which case streamline plots for Re = 46, 48, and 50 are shown. The streamlines at 
Re = 46 and 48 show symmetric vortices implying the steady nature of the flow, whereas at Re = 
50 the streamlines are asymmetric showing the unsteady nature of the flow. The time 
dependence of the flow at Re = 50 is also described in Figure 6.5 via the streamline plots 
constructed at different time steps. 

The lift coefficient history at Re = 50 calculated from QUICK scheme is shown in Figure 
6.6, which shows that the disturbance introduced into flow is sustained and resulted in periodic 
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vortex shedding from the cylinder. The streamline plots at Re = 50 for different time steps are 
presented in Figure 6.7, in which case the unsteady nature of the flow is clearly evident. The 
critical Re determined using first order upwinding scheme for momentum upwinding is 100, 
which shows more than 100% deviation from the literature value of 46 [18]. The lift coefficient 
histories at Re = 94, 96, and 100 are shown in Figure 6.8 for the first order upwinding scheme. It 
is evident from the figure that perturbation decayed at Re = 94 and 96, whereas it resulted in the 
periodic solution at Re =100. 

The comparison of and St in the range of Re from 50 to 1000 for the above- 
mentioned three-upwinding schemes is shown in Tables 6.3 and 6.4. Figure 6.9 also illustrates 
the variation in C^,. It may be seen from the figure that both second order upwinding and 
QUICK scheme predict the identical time averaged drag coefficient with a negligible variation of 
4%, whereas first order upwinding scheme predicts higher at Re < 50 and lower at Re > 

100 than predicted by the second order and QUICK schemes. The comparison of St is 
described in Figure 6.10. As seen from the figure, second order and QUICK schemes predict the 
same St for each Re and are also in good agreement with those obtained by Braza et al. [17]. On 
the other hand, first order upwinding scheme predicts lower St values than the values reported in 
the literature. 

6.2 The transition of flow past a circular cylinder confined in a plane channel from steady 
to unsteady 

As shown in Section 6. 1 , the steady flow past a circular cylinder in unbounded domain 
loses stability leading to periodic vortex shedding as Re, based on the diameter of the cylinder, is 
increased above 50. However, when the circular cylinder is confined in a plane channel, in which 
case the velocity profile is parabolic, the dynamics of the stability of the flow differ greatly 
owing to the blockage effect induced by the stationary walls of the plane channel. The effect of 
blockage on the stability of the flow may also vary depending on the situation, whether the 
cylinder is confined symmetrically between the plane walls or the cylinder is closer to one of the 
confining walls i.e. the cylinder is confined asymmetrically. The typical grids used to simulate 
flow past a cylinder confined symmetrically and asymmetrically in a plane channel are shown in 
Figure 6.1 1. The grid sensitivity analysis is done for each case to determine the optimum grid 
size. The Reynolds number for the flow past a cylinder confined in a rectangular plane channel 


46 



is defined in a way that is useful for the comparison of the results with the unbounded flow past a 
circular cylinder. 


fJ- 


( 6 . 1 ) 


where is the velocity average over the diameter of the cylinder, given as follows: 
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dy 


( 6 . 2 ) 


where dj and d 2 are the y-coordinates of the cylinder and are function of P and y as given below 




(6.3) 

(6.4) 


6.2. 1 Validation of critical Re 

The transition of flow from steady to unsteady is delayed when the circular cylinder is 
confined in a plane channel. Chen et al. [23] have determined the critical Re at which the flow 
becomes unsteady for = 0 to 0.7 for a symmetrically confined cylinder. In their study, it was 
shown that the critical Re increases with increase in f i.e. the transition is delayed with increase 
in 0. The model in the present study has been validated with their results for the critical Re and St 
for y5= 0 to 0.4 and the results are found to be in good agreement, as shown in Table 6.5. The 
comparison is also described in Figures 6.12. For determining critical Re, the flow is disturbed 
by the rotation of the cylinder as described in section 6.1.1 for unbounded flow and the ensuing 
time dependent motion is examined to check whether the disturbance is decayed or amplified 
leading to vortex shedding. Zovatto and Pedrizzetti [24] have determined the critical Re for a 
cylinder confined asymmetrically between the plane walls for 0= 0.2 and = 1 to 0.125 and the 
critical Re was found to increase with the decrease in y. We have validated our model with their 
results for critical Re and the results are in good agreement, as shown in Table 6.6 and Figure 
6.13. 

6.2.2 Variation of critical Re with 0 and y 
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As shown in section 6.2.1, the transition of the flow past a cylinder confined 
symmetrically in a plane channel is delayed with increase in p. In the present study, the variation 
of critical Re with yfor /?= 0. 1 to 0.4 is determined using time dependent numerical simulations. 
Table 6.7 presents the results for the variation of critical Re with P and y. While determining the 
critical Re, Re was increased in steps of 4 in the simulations. Therefore, the difference between 
the actual critical Re and the number determined in this study may be around 3. The effect of 
variation in y?and yon critical Re is also described in Figure 6.14. As observed from the figure, 
for a given ^the critical Re increases with decrease in y implying that the transition of the flow is 
delayed as the cylinder moves closer to one of the confining walls. It may also be observed from 
the figure that the critical Re increases steeply at y < 0.5. The increase is more pronounced at 
higher p. For example, at /?= 0.3, the critical Re increases from 120 to 276 as y decreases from 
0.5 to 0.25. This may be attributed to the fact that as the cylinder approaches the wall, the 
interaction of the boundary layer of the wall with that of the cylinder suppresses the vortex 
shedding from the cylinder resulting in the stability of the flow. Figure 6.14 also describes the 
effect of P for a constant y. The critical Re increases with increase in P for a given y, which is in 
agreement with the general trend observed that the vortex shedding is suppressed with increase 
in p. 

6.2.3 Effect of blockage ratio (P) on , St, Cl and Nuwfor a symmetrically confined cylinder 

The presence of the stationary walls significantly affects the drag force exerted on the 
cylinder by the fluid. Table 6.8 presents the effect of increase in P on for different Re for a 
cylinder confined symmetrically in the plane channel. The effect is also pictorially depicted in 
Figure 6.15. As observed from the figure, C,^ increases with increase in P from 0 (unbounded 
flow) to 0.1 for Re < 50. For Re > 50 the effect of increase in p from 0 to 0.1 on is 
negligible. It may also be observed from the figure that as the value of p increased from 0. 1 to 
0.4, C,y is increases for fixed Re. The increase inC„ with increase in ffis attributed to the fact 

that as the walls approach closer to the cylinder, the obstruction to the flow increases. Therefore, 
the velocity gradients become relatively steeper resulting in the higher drag force on the cylinder. 
However, the increase in C/, is pronounced only at low Re. The significant increase in C^j is 
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caused by the large boundary layer thickness at low Re i.e. the effect of viscous forces extends to 
lager distance from the cylinder resulting in the interaction with walls. 

Table 6.9 and Figure 6.16 presents the effect of increase in /3 on St. also illustrates the 
effect of /? on St. As seen from the figure, for all values of /?, except at 0.4, St increases with 
increase in Re, which is in agreement with the general trend. The change in St as p increases 
from 0 to 0.2 is negligible for Re < 300, whereas it decreases for Re > 400. For example, at Re = 
500 the decrease in St is 3.7% as p increases from 0 to 0.2. It may also be observed from the 
figure that as P increases from 0.2 to 0.4, St increases significantly for a fixed Re. The increase in 
St from P= 0.2 to 0.4 at Re = 400 is 31%. For /?= 0.4, St decreases with increase in Re from 200 
to 300 and then increases with Re. 

Figure 6. 1 7 describes vorticity contours at Re = 200 for various values of p. As seen from 
the figure, at /?= 0.1 vortices that are shed from the cylinder are regularly placed and the effect 
of vortices shed from the walls is negligible on those shed from the cylinder. However, at P= 0.2 
the dominance of vortices shed from the walls grows. This effect is more pronounced as p 
increases. At /? = 0.3 and 0.4, vortices shed from the walls significantly affect the vorticity 
distribution resulting in the distortion of the symmetric distribution of vortices shed from the 
cylinder. 

Figure 6. 1 8 illustrates the trajectories of Cy and Ci during one period of oscillation of the 
alternating vortices. As seen from the figure, the amplitude of oscillation of Ci decreases with 
increase in p. The double rings in the trajectories show that the time period of oscillation of Cl 
is twice that of Cy . This may also be explained from Figure 6.19 where the trajectory of Cy and 
Cl for P= ^.\ and the time histories of Co and Cl are plotted over one period. As seen from the 
figure, from point 1 to point 2 Cl increases and Co decreases. At point 2, the reverse in the trend 
of Co takes place whereas Cl continues to increase upto point 3. From point 3 to 4 Co and Cl 
both decrease. At point 4, the reverse in the trend of Co takes place whereas Cl continues to 
decrease upto point 5. The frequency of the oscillation of Co is twice that of Cl due to the 
contribution of upper and lower alternating vortices to the oscillation of Co- 

The effect of Pon Nu,,- is illustrated in Table 6.10 and Figure 6.20. It may be observed 
from the figure that the effect of increase in P from 0 to 0.3 on Aw,,, is negligible. As P increases 
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from 0.3 to 0.4, Niu increases. For example, at Re = 100 increase in Nu^xs 7.78%, 
corresponding to increase in /?from 0.3 to 0.4. 

6.2.4 Effect of gap ratio (y) on C,^ for an asymmetrically confined cylinder 

As shown in section 6.2.3, the presence of stationary walls significantly affects the drag 
force exerted by fluid on a cylinder confined symmetrically in the plane channel. Furthermore, 
when the cylinder is confined asymmetrically in the plane channel, the dynamics of flow differs 
due to the presence of stationary wall closer to the cylinder surface. Tables 6.1 1-6.14 present the 
effect of decrease in y on for f = 0.1, 0.2, 0.3, and 0.4 respectively. Figure 6.21 also 
describes the effect of y on . As observed from the figure, at = 0. 1 as decreases i.e. the 
cylinder approaches one of the confining walls for fixed Re, increases. Change in is 
negligible as decreases from 1 to 0.75, except at Re =1000 where C,y decreases by 12%. For 
fixed Re < 200, C/j increases as y decreases from 0.75 to 0.25 and the change in is 
negligible at Re > 200. For example, at Re - 20 increase inC^, is 21% as / decreases from 0.75 
to 0.25. For = 0. 1 , the effect of decrease in yon is more pronounced only at low Re and the 
effect vanishes at Re > 200. 

It may also be observed from the figure that for f = 0.2, 0.3, and 0.4 drag coefficient 
increases with decrease in y. At higher Re and higher blockage ratios increase inC^.^ as y 
decreases from 1 to 0.25 is large. For example at /?= 0.1 and Re = 100 increase in Cjr, as / 
decreases from 1 to 0.25 is 8%, whereas at 0.4 increase in is 38%. Increase in C,^ as the 
cylinder approaches closer to one of the confining walls is attributed to increase in the resistance 
to the flow. When the cylinder is present closer to the wall, the gap between the cylinder and the 
wall is filled by the boundary layers of the cylinder and wall where viscous forces are dominant. 
In this case the viscous and pressure forces exerted on the cylinder increase resulting in increase 
in the total drag coefficient. The significant increase in at low Re is due to the large 
boundary layer thickness. 

Figures 6.22 and 6.23 illustrate the effect of yon Cdp and Cdf ^or p= 0.1, 0.2, 0.3, and 
0.4 in the steady flow range. As seen from the figures, Cdp and Cdf increase with decrease in y 


50 



for fixed Re at each /?. However, increase in Cdf is higher than increase in Cdp at /? = 0.1, 
whereas for /?= 0.2, 0.3, and 0.4 increase in Cdp is higher than increase in Cdf- For example, at 
Re — 20 and ^=0.1 increase in Cdf is 85% whereas increase in Cdp is only 7% as x decreases 
from 1 to 0.25 and at Re = 20 and 0.2 increase in Cdp is 69%, whereas increase in Cdf is 
only 12%. 

6.2.5 Effect of gap ratio (y) on vortex shedding from an asymmetrically confined cylinder 

As shown in section 6.2.3, the presence of the stationary walls significantly affects St for 
a cylinder confined symmetrically in the plane channel. The asymmetric confinement of the 
cylinder in a plane channel also considerably affects St. Tables 6.15 - 6.18 present the effect of 
decrease in ^on St for various values of f. The effect of ^on St is illustrated in Figure 6.24. As 
seen from the figure, for a given fanA Re, St increases with decrease in y. At j3 =0.1, increase in 
St at Re = 200 is 14.3% as decreases from 1 to 0.25. It may also be observed from the figure 
that for a fixed yand Re, St increases with increase in fr 

Figure 6.25 illustrates the effect of /?and ^ on vortex shedding from the cylinder at Re = 
200. In the figure the instantaneous vorticty contours are presented to show the effect of 
asymmetrical confinement of the cylinder on vortex shedding. For the simulation purpose /? is 
varied from 0.1 to 0.4 and / from 1 to 0.25. As observed from the figure, at /?= 0.1, when the 
cylinder is symmetrically confined in the plane channel, the vortex street consists of positive 
(shed from the bottom) and negative (shed from the top) vortices shed from the cylinder, which 
are regularly placed. As the cylinder approaches one of the confining walls the regular structure 
of the vortices is distorted. The distortion is caused by the presence of the wall close to the 
cylinder in which case the vortices shed from the wall begins to interfere with that from the 
cylinder. It may also be observed from the figure that at y?= 0.1 and / =1 the vortex structure 
consists of double row of vortices shed from the cylinder. However, at / = 0.75, the positive 
vortices shed from the cylinder are significantly distorted due to the effect of vortices shed from 
the wall which results in the formation of vortex street with a single row of vortices. As observed 
from the figure, at / = 0.75 the vortex street consists of double row of vortices but the positive 
vortices shed from the cylinder are weaker. At y= 0.75 and 0.25, the effect is more pronounced 
with weak positive vortices shed from the bottom side of the cylinder. In this case the vortex 
street has almost become a single row of vortices. 

A— 
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The above phenomenon may also be explained by the trajectories of Ci-Cu over one 
period of the unsteady periodic state, as shown in Figure 6.26 for /?= 0.1 and Re = 200. The 
double ring in the trajectories is due to the fact that the frequency of the drag coefficient is twice 
that of the lift coefficient. Re-referring Figure 6.19, this in turn is due to the simultaneous 
contribution of the negative and positive vortices shed from the cylinder to the drag. The 
phenomenon of the formation the double ring has been explained in section 6.2.3. As observed 
from Figure 6.26, at y =1. i.e. when the cylinder is confined symmetrically in the plane channel, 
the lift coefficient oscillates symmetrically about zero with an amplitude of ±0.64. The mean lift 
coefficient in this case is zero. However, at larger values of y(0.75, 0.5 and 0.25) the oscillations 
in the values of the lift coefficient are not symmetric about zero. Therefore, the mean lift 
coefficient has a non-zero value in the asymmetrically confined case. Due to the relatively larger 
amplitude of the oscillations in the lift coefficient in the negative direction the mean lift 
coefficient is negative. It may also be observed from the figure that the amplitude of the lift 
coefficient increases in the negative direction and decreases in the positive direction as the 
cylinder moves closer to the bottom wall of the plane channel. This phenomenon is directly 
related to the shedding of the negative and positive vortices from the cylinder. The double ring in 
the trajectories of Cl-Cd gradually begins to collapse into a single ring as the cylinder moves 
closer to the wall, as seen in the figure for 0.1 in the case of y = 0.25. 

Figure 6.25 illustrates the effect of decrease in y on the vortex shedding from the 
cylinder for = 0.2 at Re = 200. As seen from the figure, when the cylinder is symmetrically 
confined in the plane channel, the structure of vortex street is well defined, with regularly placed 
positive and negative vortices. However, when compared to the symmetrically confined case at J3 
= 0.1, the vortex shedding structure is relatively irregular. It may also be observed from the 
figure that as the cylinder approaches one of the confining walls the double row of vortices 
converges to single row of vortices. At y = 0.125, i.e. when the cylinder is present very close to 
the bottom wall, the vortex shedding is completely suppressed. 

Figure 6.27 describes the suppression of vortex shedding from the cylinder in which the 
trajectories of Cl-Cd for /? = 0.2 at various gap ratios are shown. As seen from the figure the 
amplitude of Cl in the negative direction becomes relatively larger as the cylinder moves closer 
to the wall. At y = 0.125, the trajectory reduces to a point with a positive lift. In such case the lift 
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and drag coefficients are constant and are not fluctuating with time, indicating that the vortex 
shedding from the cylinder is completely suppressed. 

Re-referring Figure 6.25, the vortex shedding is gradually suppressed as the cylinder 
approaches one of the confining walls at = 0.3 and 0.4. As seen from the figure at / = 0.25 
vortex shedding from the cylinder is completely suppressed in both the cases. The 
corresponding trajectories of Cl-Cq 31/3= 0.3 and 0.4 are shown in Figures 6.28 and 6.29. These 
trajectories also confirm the suppression of vortex shedding from the cylinder. Figure 6.30 
illustrates the effect of y on Cj for various blockage ratios. As observed from the figure, for 
fixed Re, C, decreases as y decreases. The trend is reversed when the cylinder is present 
sufficiently close to the cylinder which is observed for /?= 0.2, 0.3 and 0.4. At 0.1, the 
cylinder is not present sufficiently close to the cylinder so as to results in the mean positive lift. 
The unsteady vorticity contours at Re = 200 for /3= 0.3 at various yare shown in Figure 6.31. 
The steady nature of the flow at y = 0.25 is evident from the figure. 

6.2.6 Effect of gap ratio (}j on Nuw for an asymmetrically confined cylinder 

Tables 6.19 - 6.22 present the effect of decrease in yon Num^ox fi= 0.1, 0.2, 0.3, and 
0.4, respectively. Figure 6.32 also describes the effect of yonNuw . It may be observed from the 

figure that the effect of decrease in yoniVw„- is negligible for all blockage ratios. The variation 
of local Nuy„ on the surface of the cylinder with that in the gap ratio is described in Figure 6.33. 
As observed from the figure, at each blockage ratio change in the distribution of local on the 
surface of the cylinder is negligible with decrease in y, except some distortion in the distribution 
when the cylinder is present very close to the bottom wall. The distortion is more pronounced at 
high blockage ratios. However, change in the surface averaged Nu.^ is negligible with decrease in 
y for each blockage ratio. 

6.3 The transition of flow past a circular cylinder from 2-D to 3-D in unbounded domain 

As discussed earlier, the flow past a circular cylinder in unbounded domain loses stability 
around Re = 47, leading to periodic vortex shedding. If Re is increased to 180 flow transition 
takes from 2-D to 3-D with a complex wake structure. Williamson [50] has investigated 
experimentally the discontinuities in St-Re curve in the range of Re from 49 to 178. The author 


53 



has concluded that the discontinuities in St-Re curve are caused by the three dimensional 
transition in the wake of the cylinder due to the change in oblique vortex shedding modes. The 
‘Quasi-periodic’ mode that exists for Re < 64 consists of regions of different frequencies along 
the span of the cylinder. In this mode the central region of higher frequency (fu) are sandwiched 
between two regions of lower frequency (fi). When Re is increased above 64, shedding mode 
changes to ‘Periodic’ oblique shedding mode, which consists of the regions of lower frequency 
along the entire span with a ‘chevron’- shaped pattern. In both of the shedding modes there are 
end regions of lower frequency (f^ that are significantly affected by the end conditions of the 
cylinder. It was also shown that the end regions merge together covering the whole span if the 
aspect ratio (L/D) of the cylinder is less than 28. By manipulating the end conditions of the 
cylinder, the parallel mode of vortex shedding was observed which indicates that indeed there 
exists a universal St-Re curve without discontinuities. The relation between St for the parallel and 
oblique shedding modes was given as S© = Sq/cos 9, where So and Se are Strouhal numbers for 
parallel and oblique (6) shedding modes respectively. 

The existence of two successive transitions in the transition of the flow past a cylinder to 
three-dimensionality was also observed by Williamson [51]. The first discontinuity that exists 
between Re of 170-180 is called Mode A, which has large spanwise wave lengths of 3-4 
diameters and scales to the larger physical structure in the flow i.e. the primary vortex core. The 
second discontinuity appears between Re of 230-260 is called Mode B, which has small 
wavelength of around 1 diameter and scales to the smaller physical structure of the flow i.e. the 
braid shear layer. In these modes, vortex dislocation occurs resulting in the three-dimensionality 
of the flow. It was observed by Prasad and Williamson [52] that vortex shedding control is 
possible by manipulating the end conditions of the cylinder, except in the transition region 190- 
260 which is marked by vortex dislocations. 

In the present study, we have cafried out 3-D numerical simulation to characterize the 
discontinuities in St-Re relationship and vortex dislocation. The geometry used to simulate 3-D 
unbounded flow past a circular cylinder, described schematically in Figure 3.2, is a rectangular 
brick of dimensions 28^7, 16J and 2.25d in x, y and z-directions, respectively, where d is the 
diameter of the cylinder. The cylinder is assumed to be located at \M and 18(7 from the inlet and 
outlet, respectively. The fluid enters at the plane ABCD with a uniform velocity t/jr^nd leaves at 
the plane EFGH. The typical mesh configuration used for 3-D simulation is shown in Figure 
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6.34. The grid independence test is carried out on this geometry as shown in Table 6.23. As 
observed from the table, change in C^, as the number of cells is increased from 282,000 to 
66^,000 is below 2%. Therefore, the grid with 282,000 cells was chosen as the optimum grid. 
The comparison of Cp and St values with the values reported in the literature is presented in 
Table 6.24. The maximum deviation in Cj^ values from the values reported in the literature is 

\\% at Re = 100. The values of St are matching well with the published data with a maximum 
deviation of 9% at Re= 1 00. 

A coarser grid size with 133,000 cells is also used to simulate 3-D flow past a cylinder to 
save computational time. The corresponding results obtained using this grid are presented in 
Table 6.25. In order to characterize the three dimensionality of the flow, the St-Re curves for 2-D 
flow and 3-D flow are compared in Figure 6.35. The results obtained using both coarser and finer 
grid are shown in the figure. In the case of finer grid the discontinuity in St-Re was found at Re = 
300, whereas in the case of relatively coarser grid the St-Re curve lies below the St-Re curve for 
the 2-D Flow case. The St-Re curve obtained using coarser grid is in accordance with the 
observation of Williamson [50] that the end regions of lower frequency merge together covering 
the whole span if the aspect ratio (L/D) of the cylinder is less than 28. However, when finer mesh 
is used the Mode A discontinuity was not observed, rather Mode B discontinuity was observed at 
Re = 300. 

The normalized velocity plots for Re = 100 along the span of the cylinder in z-direction at 
2D and 5D downstream locations are shown in Figure 6.36. As seen from the figure, at both the 
downstream locations the velocity profiles along the span of the cylinder are essentially the 
same. There is no variation in the velocity profiles in z-direction. Therefore, the flow at Re = 100 
is 2-D in nature. The normalized velocity profiles for Re = 200 are described in Figure 6.37. As 
observed from the figure the flow at Re = 200 is 2-D in nature. However, as observed by 
Williamson [51], the flow at Re = 200 is 3-D. The discrepancy may due to the small aspect ratio 
(L/D = 2.25) used in this study. The normalized velocity plots at Re = 300 are shown in Figure 
6.38. The 3-D nature of the flow is clearly evident from the figure, which shows considerable 
variation in profiles in z-direction. The low aspect ratio (L/D = 2.25) used in this study is found 
to be inadequate to predict the discontinuities in St-Re curve. Hence, high aspect ratio has to be 
used to predict the three-dimensional transition in the cylinder wake, which requires large 
amount of computational resources and time. 
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Rotational speed (Q) 

rad/sec 


Cl 

St 

2 

1.303 

±0.285 

0.1608 

6 

1.303 

±0.283 

0.1608 


Table 6.1: The effect of the disturbance (rotational speed) on , Ct, and St 
at Re = 100 


Upwinding Scheme 

Rec 

Rec 

Jackson (1987) 

First order 

100 


Second order 

50 

46 

QUICK 

50 

- 


Table 6.2: Comparison of critical Re observed from various upwinding schemes 


S.No 

Re 


QUICK 

First order 

Second order 

Braza et al. (1986) 

1 

40 

1.5124 

1.6086 

1.5022 

1.45 

2 

50 

1.3970 

1.5225 

1.3936 

1.40 

3 

100 

1.3059 

1.2425 

1.3027 

1.253 

4 

200 

1.3010 

1.1747 

1.2981 

1.321 

5 

500 

1.3258 

1.1551 

1.3814 


6 

1000 

1.4741 

i 

1.189 

1.4734 

1.198 


Table 6.3: Comparison of values calculated from different upvt^inding schemes 
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S.No 

Re 

St 

QUICK 

First order 

Second order 

Brazaet al. 
(1986) 



0.1182 

Steady State 

0.1212 

0.1200 

2 

100 




0.165 

3 




0.1912 

0.198 

4 

500 

0.2081 

0.1627 

0.2222 

- 

5 

1000 

0.2376 

0.1731 

0.2370 

0.22 


Table 6.4: Comparison of St values calculated from different upwinding schemes 


S.No 


Rec 

St 




Chen et al. 

This work 

Chen et al. 




(1995) 


(1995) 

B 

0 

50 

47.9 

0.121 

0.138 

2 

0.1 

50 

51.6 

0.120 

0.122 

3 

0.2 

69 

69 

0.158 

0.158 

4 

0.3 

92 

92 

0.214 

0.215 


0.4 

110 

110.1 

0.287 

0.286 


Table 6.5: Comparison of Rec and Si values at various blockage ratios 


S.No 

r 

Rec 

This work 

Zovatto and Pedrizzetti 

(2001) 

1 

1 

69 

68 

2 

0.587 

72 

75 

3 

0.375 

108 

110 

4 

0.125 

300 

300 


Table 6.6: Comparison of jRec values for various /(fi—OJ) 
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4 

40 

1.567 

1.860 

2.47 

3.643 

5 

50 

1.426 

1.681 

- 

- 

6 

60 


1.5563 

2.030 

2.914 

7 

80 

- 


1.796 

2.543 

8 

100 

1.296 

- 

1.650 

2.314 

9 

110 

- 

- 

- 

2.229 

10 

200 

1.290 

1.253 

1.421 

1.858 

11 

300 

1.321 

1.268 

1.391 

1.733 

12 

400 

1.338 

1.274 

1.384 

1.688 

13 

500 

1.340 

1.324 

1.385 

1.658 

14 

1000 

1.400 

1.368 

1.403 

1.612 


Table 6.8: Effect of blockage ratio {p} onC,^ for / = 1 (symmetrical case) 
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S.No 

Re 

St 

o 

II 

0.2 

0.3 

0.4 

1 

50 

0.1210 

Steady 

Steady 

Steady 

4 

100 

0.1608 

- 

0.2175 

Steady 

5 

110 


- 

- 

0.2870 

6 

200 

0.1914 

0.1922 

0.2260 

0.2946 

7 

300 

0.2046 

0.2028 

0.2251 

0.2721 

8 

400 

0.2100 

0.2054 

0.2279 

0.2702 

9 

500 

0.2183 

0.2140 

0.2351 

0.2725 

10 

1000 

0.2330 

0.2261 

0.2466 

0.2885 


Table 6.9: Effect of blockage ratio (fi) on St for y = l(symmetrical condition) 


S.No 

Re 

Nu„. 

o 

1! 

0.2 

0.3 

0.4 

1 

10 

1.930 

2.063 

- 

- 

2 

20 

2.526 

2.667 

2.873 

- 

3 

30 

2.969 

3.119 

- 


4 

40 

3.337 

3.497 

3.732 

4.064 

5 

50 

3.660 

3.830 

- 

- 

6 

60 


4.135 

4.385 

4.751 

7 

80 


- 

4.948 

5.345 

8 

100 

5.219 

- 

5.455 

5.880 

9 

200 

7.596 

7.592 

7.7253 

8.1813 

1 

300 

9.493 

9.326 

9.4953 

10.0998 

11 

400 

11.011 

10.916 ! 

10.9416 

11.5889 

12 

500 

12.402 

12.131 

12.1821 

12.7945 

13 

1000 

17.300 

18.980 

17.9611 

18.4201 


Table 6.10: Effect of blockage ratio {fi) on Nu^. for y = 1 (symmetrical case) 
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S.No 

Re 


r = 1 




1 






2 

20 





j) 

30 





4 

40 





5 

50 




1.734 

6 

100 



1.296 




7 

200 

1.290 



1.325 

8 

300 

1.321 



1.323 

9 

400 




1.318 

10 

500 

iiiiiiimQii|H 




11 

1000 




1.239 


Table 6.11: Effect of gap ratio (j) onC„ for 0.1 


S.No 

Re 


II 



0.125 

1 

10 

4.175 

4.197 

4.642 

5.549 

2 

20 

2.671 

2.686 

3.149 

3.918 

o 

J 

30 

2.142 

2.154 

2.556 

3.290 

4 

40 

1.860 

1.871 

2.222 

2.940 

5 

50 

1.681 

1.691 

2.005 

2.714 

6 

60 

1.556 

1.564 

1.852 

2.550 

7 

200 

1.253 

1.2573 

1.333 

1.750 

8 

300 

1.268 

1.2676 

1.318 ^ 

1.586 

9 

400 

1.274 

1.2856 

1.380 

1.609 

10 

500 

1.324 

1.300 

1.380 

1.606 

1 1 

1000 

1.368 

1.400 

1.4305 

1.897 


Table 6.12: Effect of gap ratio WonC/, for ^=0.2 
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S.No 

Re 

c. 

r = \ 

0J5 

0.5 

0.25 

1 

20 

3.759 

3.853 

4.142 

- 

2 

40 

2.476 

2.540 

2.820 

2.726 

o 

J 

60 

2.030 

2.079 

2.308 

2.330 

4 

80 

1.796 

1.836 

2.033 

2.098 

5 

100 

1.650 

1.682 

1.859 

1.940 

6 

200 

1.421 

1.438 

1.547 

1.516 

7 

300 

1.391 

1.404 

1.484 

1.399 

8 

400 

1.384 

1.396 

1.452 

1.352 

9 

500 

1.385 

1.393 

1.435 

1.323 

10 

1000 

1.403 

1.397 

1.395 

1.266 


Table 6.13: Effect of gap ratio (jj on for ;9= 0.3 


S.No 

Re 


r = 1 

0.75 

0.5 

0.25 

1 

40 

3.643 

3.761 

4.1253 

4.390 

2 

60 

2.914 

3.004 

3.3468 

3.779 

3 

80 

2.543 

2.618 

2.9173 

3.434 

4 

100 

2.314 

2.381 

2.6482 

3.204 

5 

110 

2.229 

2.292 

- 

- 

6 

120 

- 

2.218 

2.463 

3.026 

7 , 

200 

1.858 

1.904 

2.098 

2.5716 

8 

300 

1.733 

1.748 

1.9241 

2.3044 

9 

400 

1.688 

1.683 

1.8363 

2.1714 

10 

500 

1.658 

1.658 

! .7974 

2.0868 

1 1 

1000 

1.612 

1.620 

1.7356 

1.8986 


Table 6.14: Effect of gap ratio onC^j for 0.4 
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S.No 

Re 

St 

y = 1 

0.75 

0.5 

0.25 

1 

50 

0.1210 

Steady 

Steady 

Steady 

2 

55 

- 

0.1273 

Steady 

Steady 

3 

100 

0.1608 

0.1610 

0.1660 

0.1973 

4 

200 

0.1914 

0.1883 

0.1930 

0.2188 

5 

300 

0.2046 

0.2023 

0.2055 

0.2233 

6 

400 

0.2100 

0.2120 

0.2108 

0.2266 

7 

500 

0.2183 

0.2170 

0.2170 

0.2287 

8 

1000 

0.2330 

0.2377 

0.2380 

0.2440 


Table 6.15: Effect of gap ratio (/) on St for >9= 0.1 


S.No 

Re 

St 

r = 1 

0.875 

0.375 

0.125 

1 

69 

0.1589 

Steady 

Steady 

Steady 

2 

72 


0.1614 

Steady 

Steady 

3 

108 

- 

- 

0.2142 

Steady 

4 

200 

0.1922 

0.1942 

0.2346 

Steady 

5 

300 

0.2028 

0.2051 

0.2349 

0.3063 

6 

400 

0.2054 

0.2114 

0.2361 

0.2938 

7 

500 

0.2140 

0.2148 

0.2377 

0.2608 

8 

1000 

0.2261 

0.2268 

0.2300 

0.2843 


Table 6.16: Effect of gap ratio (j) on Si for 0.2 
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S.No 

Re 

St 

r = i 

0.75 

0.5 

0.25 

1 

100 

0.2175 

0.2205 

Steady 

Steady 

2 

200 

0.2260 

0.2341 

0.2563 

Steady 

3 

300 

0.2251 

0.2320 

0.2527 

0.26388 

4 

400 

0.2279 

0.2330 

0.2488 

0.27037 

5 

500 

,0.2351 

0.2358 

0.2483 

0.26979 

6 

1000 

0.2466 

0.2460 

0.2555 

0.27097 


Table 6.17: Effect of gap ratio (y) on St for 0.3 


S.No 

Re 

St 

r = 1 

0.75 

0.5 

1 

no 

0.287 

Steady 

Steady 

2 

120 

- 

0.2852 

Steady 

4 

200 

0.2946 

0.3006 

0.3089 

5 

300 

0.2721 

0.2791 

0.3113 

6 

400 

0.2702 

0.2684 

0.3058 

7 

500 

0.2725 

0.2687 

0.3046 

8 

1000 

0.2885 

0.2790 

0.2926 


Table 6.18: Effect of gap ratio (]) on St for 0.4 
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S.No 

Re 

Nuw 

r = \ 

0.75 

0.5 

0.25 

1 

10 

1.930 

1.9281 

1.9421 

1.9039 

2 

20 

2.526 

2.525 

2'.5501 

2.629 

3 

30 

2.969 

2.9686 

3.0014 

3.1107 

4 

40 

J.JJ / 

3.3365 

3.3756 

3.5038 

5 

50 

3.660 

3.6598 

3.70433 

3.8478 

6 

100 

5.219 

5.2077 

5.2271 

5.2927 

7 

200 

7.596 

7.5405 

7.6368 

j 

7.7468 

8 

300 

9.493 

9.3808 

9.3426 

9.4156 

9 

400 

11.011 

10.935 

10.904 

11.105 

10 

500 

12.402 

12.1702 

12.212 

12.6958 

11 

1000 

17.300 

18.291 

18.005 

18.6624 


Table 6.19: Effect of gap ratio (j) on //w« for J3- 0.1 


S.No 

Re 

Nuw 

r = 1 

0.875 

0.375 

0.125 

1 

10 

2.063 

2.064 

2.093 

2.755 

2 

20 

2.667 

2.670 

2.706 

3.126 

3 

30 

3.119 

3.123 

3.205 

3.435 

4 

40 

3.497 

3.501 

3.609 

3.719 

5 

50 

3.830 

3.835 

3.956 

3.994 

6 

60 

4.135 

4.140 

4.269 

4.268 

8 

200 

7.592 

7.601 

7.615 

7.520 

9 

300 

9.326 

9.354 

9.452 

9.673 

10 

400 

10.916 

10.839 

11.237 

11.415 

11 

500 

12.131 

12.088 

12.647 

13.209 

12 

1000 

18.980 

17.927 

19.336 

20.540 


Table 6.20: Effect of gap ratio ()) on for 0.2 
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S.No 

Re 

Nu „■ I 

7 = 1 

0.75 

0.5 

0.25 

1 

20 

2.873 

2.880 

2.861 

- 

2 

40 

3.7385 

3.747 

3.786 

3.610 

o 

J 

60 

4.385 

4.404 

4.461 

4.204 

4 

80 

4.948 

4.969 

5.033 

4.765 

5 

100 

5.455 

5.478 

5.546 

5.270 

6 

200 

7.725 

7.759 

7.759 

7.281 

7 

300 

9.495 

9.609 

9.595 

9.000 

8 

400 

10.941 

11.172 

11.085 

10.561 

9 

500 

12.182 

12.245 

12.413 

11.922 

10 

1000 

17.961 

17.823 

18.207 

16.560 


Table 6.21: Effect of gap ratio on Nu^ for 0.3 


S.No 

Re 

Nu^. 

r = 1 

0.75 

0.5 

0.25 

1 

40 

4.064 

4.069 

4.035 

4.160 

2 

60 

4.751 

4.762 

4.773 

4.697 

o 

J 

80 

5.345 

5.358 

5.387 

5.220 

4 

100 

5.880 

5.895 

5.937 

5.732 

5 

110 

- 

6.145 

- 

- 

6 

120 

- 

6.388 

6.441 

6.223 

7 

200 

8.181 

8.173 

8.263 

7.993 

8 

300 

10.099" 

10.076 

10.236 

9.716 

9 

400 

11.588 

11.730 

11.711 

11.156 

10 

500 

12.794 

12.903 

13.181 

12.492 

1 1 

1000 

18.420 

18.217 

18.352 

17.751 


Table 6.22: Effect of gap ratio {)) on Nuw for fi- 0.4 


65 





Re 


St 

c. 

Persillon 
and Braza 
(1998) 

St 

Persillon 
and Braza 
(1998) 

GRIDl 

GR1D2 

GRIDS 

GRIDl 

GRID2 

GRIDS 

100 

1.386 

1.391 

1.394 

0.168 

0.167 

0.168 

1.240 

0.164 

200 

1.371 

1.390 

1.369 

0.195 

0.193 

0.197 

1.306 

0.181 

500 

1.259 

1.264 

1.284 

6.199 

0.201 

1 

0.201 




Table 6.23: Grid independence test for 3-D flow past a cylinder in unbounded 
domain 


GRIDl : 282,000 cells with 96 intervals on the surface of the cylinder. 
GRID2: 461,000 cells with 120 intervals on the surface of the cylinder. 
GRIDS: 663,000 cells with 120 intervals on the surface of the cylinder. 


S.No 

Re 


St 

This 

work 

Persillon 
and Braza 
(1998) 

Zhang 

et.al. 

(1998) 

This 

work 

Persillon 
and Braza 
(1998) 

Zhang 

et.al. 

(1998) 

1 

100 

1.386 

1.240 

1.32 

0.168 

0.164 

0.154 

2 

200 

1.371 

1.306 

1.28 

0.195 

0.181 

0.190 

3 

220 

1.373 

1.313 

- 

0.199 

0.184 


4 

300 

1.366 

1.366 

- 

0.199 

0.206 

“ 

5 

400 

1.375 

- 



1 

- 

6 

500 

1.259 

- 


0.199 

- 



Table 6.24: and St for 3-D flow past a cylinder using GRID2 
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S.No 

Re 


St 

This work 

Persillon and 
Braza (1998) 

This work 

Persillon and 
Braza (1998) 

1 

50 

1.595 


Steady 

- 

2 

100 

1.336 

1.240 

0.147 

0.164 

o 

150 

1.266 

- 

0.156 

0.175 

4 

180 

1.239 


0.159 

0.178 

5 

200 

1.241 

1.306 

0.165 

0.181 

6 

250 

1.180 

1.343 

0.166 

0.186 

7 

300 

1.123 

1.366 

0.170 



Tabic 6.25: and St for 3-D flow past a cylinder obtained using coarser grid 


(133,000 cells) 
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Figure 6.1 : Typical mesh configuration used to simulate flow past a cylinder in unbounded domain (R/r = 140) 





100 150 200 250 


t = tUa/D 


Figure 6.3: Co and Cl 














Figure 6.5: Streamline plots at Re = 46, 48 and 50 for second order upwinding scheme 
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Lift Convergence History (Time-1. 3222 e+01) 


Aug 13, 2004 

FLUENT 6.0 (2d. dp. segregate^j, lam, unsteady) 


Figure 6.6; Cl vs time tor Jte = 50 using QUICK scheme 
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Figure 6.9: Cd vs Re calculated from various upwinding schemes 









Figure 6.11: Typical mesh configurations used to simulate flow past a cylinder confined in plane channel 








Figure 6.12: Comparison of Rec and St for various p and y = 1 (symmetrical) 
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Figure 6.31: Unsteady vorticity contours at Re = 200 and p = 0.3 
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CHAPTER? 

CONCLUSIONS AND SCOPE FOR THE FUTURE STUDY 

7.1 Conclusions 

The 2-D unsteady momentum and heat transfer from an asymmetrically confined 
circular cylinder in a plane channel was studied numerically. In addition, 3-D unsteady 
flow past a circular cylinder in unbounded domain was studied to characterize three- 
dimensional transition in the wake of the circular cylinder. The equations of continuity, 
momentum and energy equations were solved using ELUENT, a flow modeling software 
that uses a finite volume based numerical method. First, the numerical simulations were 
carried out for flow past a circular cylinder in unbounded domain in steady and unsteady 
flow range and the simulation results were compared with the published work. In 
particular, the transition of the flow from steady to unsteady is studied with the aid of 
various momentum upwinding schemes available in FLUENT. The global parameters 

like; Co,Nuw and St obtained using various upwinding schemes were compared. 
Simulations were also carried out for the flow past a symmetrically confined cylinder in a 
plane chaimel. In the case of cylinder confined symmetrically in a plane channel the 
transition of flow from steady to unsteady was studied. The critical Re at the flow 
transition takes place was determined for various blockage ratios (0.1 <fi< 0.4). The 
results were compared with the published work and were found to be in good agreement. 

The simulations were then extended to study the momentum and heat transfer 
from an asymmetrically confined circular cylinder in a plane channel in the range of Re 
from 10 to 1000 and 0.1 < ^< 0.4. In this case, the transition of flow from steady to 

unsteady is determined as a function of and /. The effect of and yon and St 

is studied. It is observed that for a symmetrically confined cylinder the critical Re 
increased with increase in p. For fixed Re, and St increased with increase in fi. The 

effect of P on Nuw was found to be negligible for p < 0.3. For a cylinder confined 
asymmetrically in a plane channel, the critical Re at which the flow transition takes place 
was found to increase with decrease in gap ratio (j), for fixed P, i.e. the vortex shedding 
from the cylinder is suppressed as the cylinder approaches one of the confining walls. 
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The critical Re was also found to increase with increase in piov fixed y- For fixed Re, 

was found to increase with decrease in y i.e. the drag force exerted on the cylinder 
increases as the cylinder approaches closer to the wall. The increase in drag coefficient 
was more pronounced at low Re. St was found to increase with increase in Re for fixed p 
and y It was also found that, for fixed Re St increases with decrease in y. When the 
cylinder is confined symmetrically in a plane channel, the stmcture of vortices shed firom 
the cylinder is well defined with double row of vortices placed regularly. The double row 
becomes a single row of vortices that are shed fi'om the top of the cylinder as the vortex 
shedding firom the bottom of the cylinder is suppressed due to the interference of the 
vortices shed form the bottom wall. The vortex shedding firom the cylinder is completely 
suppressed when the cylinder is present very close to the cylinder. For fixed Re, the 
amplitude of the oscillations in the lift coefficient decreases as p increases. For fixed Re 
and p, the oscillation in the lift coefficient increases in the negative direction as y 
decreases. The trend is reversed when the cylinder is present sufficiently close to the 
cylinder i.e. is positive when the cylinder is present very close to the cylinder. The 

effect of decrease in y on Nuw was found to be negligible for all blockage ratios. 
However, some distortion was found in the distribution of local Nu^ on the surface of the 
cylinder as the cylinder moves closer to the wall. 

7.1 Scope for the future study 

(i) The present study was carried out assuming laminar flow conditions. However, 
turbulent flow may also be simulated for high range of Re. 

(ii) In the case of flow past an asymmetrically confined cylinder, constant wall 
temperature boundary condition was assumed. Constant wall flux boundary may be used 
to find the effect of asymmetry on the haet transfer from the cylinder. 

(iii) In the present work, an incompressible Newtonian fluid was considered. Therefore, 
the simulations may be carried out using non-Newtonian fluid as well. 

(iv) In the present work, cross flow situation was considered. This may be extended to 
parallel flow past a cylinder in unbounded and bounded domains. 
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APPENDIX 


This section presents the procedural steps needed to be followed in order to 
solve a particular type of flow problem using the flow modeling software, FLUENT. 

1. In a UNIX environment, one has to type gambit, followed by a session name (for 
example modell) in the console panel. Then GAMBIT is opened and modell is started. 

2. The desired geometry (2D or 3D) is created using face/volume creation panel. 

3. The flow domain consisting of that particular geometry is now required to be meshed. 
A 2D geometry can be meshed using triangular or quadrilateral elements. However, for 
a 3D geometry, meshing can be achieved by using hexahedral, wedge-shaped or 
tetrahedral elements. It may be pointed out that each element in GAMBIT is associated 
with a set of type options. Map, submap, and pave are a few examples of available type 
options. For example, meshing with triangular elements can be achieved by the type pave, 
which allows for unstructured meshing throughout the domain. 

4. Once the meshing is completed, proper boundaries are to be specified from the 
boundaries specification panel for each edge/face; and continuum type (fluid or solid) is 
also specified for face (2D) or volume (3D). 

5. Then the existing file is saved as a mesh file using File/Save... menu item. 

6. This saved file is then exported to FLUENT using File/Export. . .option. 

7. Afterwards FLUENT is opened by typing fluent in the console panel. In the fluent 
panel proper solver version (2D or 3D etc.) has to be specified in the command line. 
Then from the File/Read. . . menu the saved mesh file is being read. 

8. At this point it is necessary to check the grid from Grid/Check... option. An error 
occurs for an improper or incomplete meshing. In such a situation the grid check fails. 

9. Fluent provides three different solver formulations: (i) segregated, (ii) coupled 
implicit, and (iii) coupled explicit. From Define/Model/Solvfeir... panel a suitable solver 
formulation is selected. From the same panel steady or utifeteady condition is selected 
considering the time-dependence of the problem to be solved. 

10. The Defme/Materials... panel is then opened and propertied of the chosen material are 
inserted manually or selected from the materials drop-down list. 



11. Afterwards the reference values panel is opened from the Define/Reference 
Values... item and the reference pressure location is selected. If gravity is important in 
the problem, the magnitude and direction of the gravitational acceleration are being 

inserted into this panel. 

12. The specification of boundary conditions remains a very important aspect in the 
simulation studies. From the Define/Boimdary Conditions... panel, proper boundary types 
are selected and corresponding boundary conditions are being set. 

13. Now, from the Solve/Controls/Solutions... panel, proper discretization schemes for 
pressure, pressure-velocity coupling and momentum are selected and under- relaxation 
factors are adjusted. 

14. At this point the convergence criteria for the residuals of momentum and velocity 
components (for isothermal and laminar flow situations) are set using the 
Solve/Monitors/Residual... panel. The plot option is also activated in order for the 
convergence history to be visible while iterations are going on. 

15. Then it is necessary to store all the input parameters in a case file (for example 
model 1 .cas) in order for the information to be available in a next FLUENT session. It is 
done from the FileAVrite/Case. . . menu item. 

16. The solution initialization is another important step that is to be completed before 
start of the iterations. Tlus initialization is achieved from the Solve/Initialization... panel 
where one can retain default values for velocity and pressure (which are zero) or select a 
zone from the Compute From... drop down list; and can specify the values of the 
pertinent variables. 

17. From the Solve/Iterate... panel the computations are now started by specifying the 
number of iterations and clicking on Iterate. However, if unsteady problem is to be solved 
then in the iterate panel another two important parameters need to be specified also. 
These are the time step size and the number of time steps. Here, instead of mentioning 
the number of iterations, one has to input the maximum number of iterations per time 
step. 

18. When the computations are being carried out, the values of the residuals of 
momentum and velocity components are .observed to decrease if all the input 
parameters are set properly and a realistic model being chosen. The iterations are 



observed to converge after a certain time when the convergence criteria are being 
satisfied. 

19. After convergence is achieved, it is necessary to save the information for further 
processing. This is achieved by FileAVrite/Data... option and the data file named 
modell.dat is being created. 

20. The converged solution enables one to obtain the values of many important 
parameters and produce graphical representations of various functions. 

21. As a representative case one might compute the value of the total drag coefficient, 
C», which is an important engineering parameter in the study of flow past submerged 
objects. This can be achieved firom the Report/Reference Values... panel where one has to 
specify the velocity, density and the projected area for the computation of Co. Now the 
Report/Forces. . .panel is to be opened and the proper zones are to be selected. Afterwards 
clicking on the Compute... button will give the drag force components and the drag 
coefficient. 

22. It is also important to study the vectors of velocity, contours of stream function, etc., 
to properly understand the flow behavior. This can be achieved from the 
Display/ Velocity Vectors.,. or Display/Contours... panel, respectively. The surface 
vorticity is another important parameter, which helps to compute the angle of separation 
(the angle at which the flow separation occurs from the surface of any submerged 
object). This can be achieved from the XY Plot... option where the surface vorticity may 
be plotted as a function of curve length. 

The above section briefly discussed several steps which are to be followed 
to simulate a standard problem of fluid dynamics using a very sophisticated flow 
modeling software, FLUENT. The readers are requested to consult the manuals for a 
better understanding of different features available in this package 



